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PKEFAOE 

I T is now over thirty years since Heaviside's operational methods of 
solving the differential equations of physics were first published, but 
hitherto they have received very little attention from mathematical 
physicists in general. The chief reason for this lies, I think, in the lack 
of a connected account of the methods. Heaviside's own work is not 
systematically arranged, and in places its meaning is not very clear. 
Bromwich's discussion of his method by means of the theory of functions 
of a complex variable established its validity ; and as a matter of practical 
convenience there can be little doubt that the operational method is far 
the best for dealing with the class of problems concerned. It is often said 
that it will solve no problem that cannot be solved otherwise. Whether 
this is true would be difficult to say; but it is certain that in a very 
large class of cases the operational method will give the answer in a 
page when ordinary methods take five pages, and also that it gives the 
correct answer when the ordinary methods, through human fallibility, 
-•re liable to give a wrong one. In particular, when we discuss the small 
jscillations of a dynamical system with n degrees of freedom by the 
method of normal coordinates,, we. obtain a determinantal equation of 
the nth degree to give the speeds of the normal modes. To find the ratios 
of the amplitudes we must then complete the solution for each mode. If 
we want the actual motion due to a .given initial disturbance we must 
solve a further family of 2 n simultaneous equations, unless special 
simplifying circumstances are present. In the operational method a 
formal operational solution is abtained with the same amount of trouble 
as is needed to give the period- equation in the ordinary method, and 
from this the complete solution is obtainable at once by a general rule 
of interpretation. For continuous systems the advantage of the opera- 
tional method is even greater, for it gives both periods and amplitudes 
easily in problems where the amplitudes cannot be found by the ordinary 
method without a knowledge of some theorem of expansion in normal 
functions analogous to Fourier's theorem. Heat conduction is also 
especially conveniently treated by operational methods. 

Since Bromwich’s discussion it has often been said that the operational 
method is only a shorthand way of writing contour integrals. It may be ; 
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but at least one may reply that a shorthand that avoids the necessi' 
of writing jr- I dXin every line of the work is worth while. Connect* 

-oTTt Jc- too 

with the saving of writing, and perhaps largely because of it, is the fe 
that the operational mode of attack seems much the more natural win 
one has any familiarity with it. After all, the use of contour integra 
in this connection was introduced by Bromwich, who has repeatedly c? 
dared that the direct operational method of solution is the better 
the two. 

My own reason for wilting the present work is mainly that I have 
found Heaviside’s methods useful in papers already published, and shah 
probably do so again soon, and I think that an accessible account of then, 
may be equally useful to others. In one respect I must offer an apolog\ 
to the reader. Heaviside developed his methods mostly in relation to tl i 
theory of electromagnetic waves. Having myself no qualification to wri • 
about electromagnetic waves I have refrained from doing so ; but as ■ ,r 
operators occurring in the theory of these waves are mostly of ty i 
treated here I think the loss will not be serious. It can in any case r, 
remedied by reading Heaviside’s works or some of the papers in the h 1 
at the end of this tract. I 

A chapter on dispersion has been included. The operational solution 
can he translated instantly into a complex integral adapted for evalua- 
tion by the method of steepest descents ; a short account of the latter 
method has also been given, because it is not at present very accessible, 
and is often incorrectly believed to be more difficult than the method of 
stationary phase. Two cases where the Kelvin first approximation to the 
wave form breaks down are also discussed. <>* 

My indebtedness to the writings of Dr Bromwich is evident from t ? 
references in the text. In addition, the problems of 4.4 and 4.5 are tak o * 
directly from his lecture notes, and several others are included larg ; 
as a result of conversations with him. \ 

My thanks are also due to the staff of the Cambridge Uni vers' ^ 
Press for their care and consideration during publication. . 

HAROLD JEFFREY 

St John’s College, 

Cambridge, 

1927 July 19. 


PREFACE TO SECOND EDITION 


r iiE chief change in this edition is the return to Heaviside’s notation 
of p , which I had replaced by <r, for reasons that appeared to me 
k equate at the time. But my course was approved by nobody else, and 
i ave not thought it worth while to maintain it. In places I have 
I- dified the discussion, and more stress has been laid on the notion of 
°* ite integration as the fundamental operator. 

^ ,e chapter on Bessel Functions has been rewritten, and a discussion 
£]., 3 transmission of a simple type of telegraphic signal along a cable 
t] been added. 

l _hanks are due to several friends for pointing out corrections to the 
:st edition, especially Dr S. Goldstein and Professor H. S. Carslaw. 

Dr Goldstein has in addition read the whole of the present edition in 
proof and suggested many improvements. 

My previous thanks to the staff of the Cambridge University Press 
'mst be renewed. 
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OHAPTEE I 

FUNDAMENTAL NOTIONS 


1.1. Let us consider the linear differential equation of the first order 


dy 

dx 


= By + 8, 


Cl) 


where R and S are known functions of x i bounded and integrable when 
Suppose further that y=y Q when x = 0. Let Q denote the 
operation of integrating with regard to x from 0 to x, so that 

Qy= ( ydx. (2) 

Jo 

Perform the operation Q on both sides of the equation (1). Then we 
find 

y-yo =Q(By + S), (3) 

and the right side vanishes with x. This can be rewritten in the forms 

(l-QB)y=zy 0 + QS, (4) 

y=y 0 +QS+QRy. (5) 

These are both equivalent to the original differential equation (1) to- 
gether with the given terminal condition. When we write QBy we mean, 
of course, that R is to be multiplied into y and the product integrated 
with regard to x from 0 to x. But the whole expression for y may be 
substituted in the last term of (5), giving in succession 

y-y o+ Q& + QR (2/0 + Q8+ QBy ) 

=y 0 + QS+ QR(3, 9 + QS) + QRQRfa+QS+ QRy) 

= (yo + QS) + QR (yo + Q8) + QRQR (y. + QJ3) 

+ QRQRQR(y*+QS)+ (6) 

on repeating the substitution indefinitely. We have to show that the 
infinite series (6) converges and that it is the correct solution. In 
evaluating each term R is supposed to be multiplied into the whole 
expression after it, and Q to operate on the whole expression after it* 
Suppose that within the range considered ^ . 

| R\*A> | V (7) 

where A and B are finite. Then the absolute value of the second term 


j 


1 
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is less than ABx, that of the third than A*Bx 2 /2 1, that of the fourth 
than A 3 Ba?/3U while the general term is less than A n Bx n ln\. Every 
term in the series is therefore less than the corresponding term in the 
series for B exp Ax , and the series is therefore convergent. It therefore 
represents a definite function, and on substituting it in (5) we see that 
the equation is satisfied for all values of x ; also the solution reduces to 
go when x~0 y as it should. Thus (6) is the correct solution, 

The solution can also be written 

y = (1 + QR + QRQR + QRQRQR + . . •) (^o + QS) (8) 

The operator between the first pair of brackets is the binomial expansion 
of (1 - QR)~\ carried out as if QR was merely a number. Since y Q + Q8 
is a determinate function, we can write the solution in the form 

; (») 

provided that y 0 + QS is evaluated first and that the operator (1 - QR) -1 
is expanded by the binomial theorem before interpretation. In fact (9) 
is merely a shorthand rule for writing (8). But on returning to (4) we 
see that (9) is also the solution of (4) carried out as if 1 - QR was a mere 
number. 

It is evident from (8) that the values of S for negative values of x do 
not affect the solution provided y Q is kept the same. Suppose then that 
8 was zero for all negative values of x } and that y was zero when x = 0. 
The solution would be 

y — (1 + QR + QRQR + ...) Q8, (10) 

and this solution would be unaltered if the lower limits of all the 
integrals were replaced by - 00 . But if now we add to 8 a, constant y 0 /£ 
for all values of x between 0 and f , QS will be increased by y 0 for all 
values of x greater than £, and (10) will be converted into (8). If then £ 
tends to zero, y 0 remaining the same, will tend to y 0 > and we recover 
the solution (8) with the original initial condition. The physical interest 
of this result is that it corresponds to our notions of causality. Suppose 
that the independent variable x is the time, and that y represents the 
departure of some variable from its equilibrium value; then R represents 
a property of the system and 8 an external disturbing influence. If the 
system was originally in its equilibrium state, the form (10) exhibits 
the disturbance produced by the external influence after it enters. If it 
was undisturbed up to time zero, the part of (8) depending on 8 repre- 
sents the effect of the finite disturbances acting at subsequent times, 
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while the part depending on y Q represents the effect of the impulsive 
disturbance at time 0 required to change y suddenly from zero to y$. If 
we like we can separate the solution into two parts 

(l-QB)-*y 9 and (1-QBy'QS, 

and say that the first represents the effect of the initial conditions and 
the second that of the subsequent disturbances. 


1.11* The method just given can be extended easily to cover many 
equations of the second and higher orders. Thus if our equation is 

w 

with y=y* and when x= 0, we find by integration 

%-Jh + .....(2) 

y - 2Jo- QQ(Ry + S) (3) 

This leads to the solution in a series 

y = (1 + Q 2 R + Q^RQ^R + • - ■) ( ?/« + xyi + ( 4 ) 

where (#) = [ f f(t ) d£dt (5) 

Jo Jo 

1.12. As a special example, consider the equation 

<*> 

with y = 1 when # = 0. Then carrying out the process of 1.1 (6) we get 
y = (l + aQ + cl 2 Q 2 +• . 1 (2) 

= 1 + 0# + —y + (3) 

the ordinary expansion of exp ax. Or take the equation 

®s(*£) +fl * S3Q ' (4) 

with y~l and dy/dx = 0 when x = 0. We can infer 

( 5 ) 

« 

the ordinary expansion of J Q (ar). 
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1.2. The foregoing method is due originally to J. Caqu£*'; it is a 
valuable practical method of obtaining numerical solutions of linear 
differential equations. Its extension to equations of any order, or to 
families of simultaneous equations of the first order, is more difficult 
unless special simplifications enter. But if the equations have constant 
coefficients, which is an extremely common case in physical applications, 
a considerable development is possible. This arises from the fact that 
the operator Q obeys the fundamental laws of algebra. Thus if a is a 


constant and u and v are known functions of x } 

Q(m) = a (Qu), ( 1 ) 

Q(u + v)=Qu+Qv, ( 2 ) 

Q m Q n u=Q n Q m u= Q m+n u ( 3 ) 

Consequently Q behaves in algebraic transformations just like a number. 
If for instance we have two operators f(Q) and g (Q) both expressible 
as sums of integral powers of Q, thus 

f(Q)=^o + a 1 Q + a a Q 2 ^a 3 Q 3 +... 9 ( 4 ) 

9(Q)~^q + biQ + b 2 Q‘ 2 + b 3 Q? + ( 5 ) 


where the at s and 6’s are constants, let us for a moment replace Q by a 
number z small enough to make the series converge absolutely. Form 


the product series 

f(z)g(z) = (a 0 + ctr L z + a 2 2 ?+ ...)(&<, + b 1 z + b 3 ^+ ...) 

= c 0 + c 1 z + c a z 2 + (6) 

say. Consider first the case where the series are both polynomials. Then 
if 8 is an integrable function of x 

XQ)g(Q)S<<h+c l Q+c*Q+...)s. (7) 

For the left side means 

(#o + (x>i Q + ct/2 Q 2 (bo 8 + bi Q8 + b% Q 2 8 + . . . ) 

= a 0 (b 0 8+ h QS+ b a Q 2 8+...) + a 1 Q (b 0 8+ b x Q8 + ...) 

+ #2 Q a (boS +...)+■ (8) 


* Liouville’s Journal , (2), 9 (1864), 185-222. Further developments are given by 
Fuchs, Ann. d. Matem (2), 4 (1870), 36-49 ; Peano, Math. Ann., 32 (1888), 450-466; 
H. F. Baker, Proc . Lond. Math. Soc (1), 34 (1902), 347-360; (1), 35 (1902), 333- 
378; (2), 2 (1905), 293-296; Phil . Trans., A, 216 (1916), 129-186. Caqu6 considers 
only a single differential equation, but notices that the operator in the result is 
in the form of a binomial expansion. For other operational methods based on these 
principles, but applicable to equations of order higher than the first, or to families 
of equations of the first order, the above papers of Prof. Baker may be consulted. 
Physical applications are given by W. L. Oowley and H. Levy, Phil . Mag., 41 
(1921), 584-607; Jeffreys, Proc . Lond. Math. Soc., (2), 23 (1924), 454 and 465; 
M.N.P.A.S. , Geoph. Suppl. 1 (1926), 380-383. 




B’tnSTDABIBSrTAL NOTIOITS 


Then using (3) we can collect the terms involving the same power of Q, 
and by (1) and (2) we can rearrange the series and obtain 

a„b 0 S + (« 0 &i + «A) QS + (fflo&a + a 1 b l + <* 2 & 0 ) + •••> •••(9) 

which is by definition the same as 

(co+CiQ + CjiQ 2 + ...)& (10) 

When the series are infinite, their convergence may be established 
easily. If the series for f(z) converges absolutely when \z\ ^r, then a 
number M must exist such that a n r n ^ M for all positive integral values 
of n. Thus 

( 11 ) 

and also if for all values of x, | 8\ ^ G, where G is a constant, 


Q n 8^G* 


Hence if 


f(Q)S=%a n Q»S, (13) 


each term is less than the corresponding term of the series 

S f®\ n G 


3 Ml- 
*~0 \r 


which converges for all values of x however large. So long as f{z) is 
expansible about the origin in a convergent power series, however 
small its radius of convergence may be, the expression j(Q)8 will be 
the sum of an absolutely convergent series however great x may be. If 
f(z) and g (z) are both expansible within some circle, their product 
series will also converge within this circle and the expression/(Q) g ( Q ) 8 
will be an absolutely convergent series however great x may be. 
Provided with this result we can now easily extend the result (7) to the 
case where/ (z) and g(z ) are infinite series, by methods analogous to 
those used to justify the multiplication of two absolutely convergent 
power series. 

1.3. We can now extend these methods to the solution of a family 
of n simultaneous differential equations of the first order with constant 
coefficients. Suppose that the equations are 

011 X \ 4 " <?12 x % 4 * . * . + &IJI X n — Si 1 

0qi Xi "h OpiXn + ••• 4* 02n X n = Sq I / 1 *\ 


0 qi X \ + 0 ‘ 2‘j X % 4 * 


4 " 4 - ... 4 - 0 ^ X n — 8 % 
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where the afa are the dependent variables, the independent variable is 

d 

t, and e rs denotes a r3 ^ + b r8j where a T8 and b r8 are constants; the are 
known functions of t. We do not assume that 


0h r8 — (l ST) b rs — bgr, 

but we do assume that the determinant formed by the a r8 is not zero. 
The terminal conditions are that when t = 0, x l — and so on, where the 
u’s are known constants. 

The typical equation (1) can be written shortly 

^t 8 0 r8 fl/ 8 ^ $ r (j $ == 1 , 2 , ( 2 ) 

or, using the convention familiar in tensor calculus, that a repeated suffix 
such as s is to be given all possible values and the results added, simply 

6r 8 &s ~ * (3) 

First perform the operation of definite integration with regard to t 
from 0 to t on both sides of this equation. We have 


Qe r , x, = 


dt 


= a r& + brs 


— fra *&8 ~ m Mrs % a j (4) 

where/™ denotes the operator a r8 + b r8 Q. Then the original equations 
and the initial conditions together are equivalent to the n equations 


fra 9D a —' f Vr^ r (5) 

where v r -a ra u 8 (6) 


Now let B denote the operational determinant formed by the f ra , 
namely 


II/JI- 


fn fi a fm 

/21 ^22 ••• fn 

fnl fn 2 • • ■ fnn 


•(7) 


If this determinant is expanded by the ordinary rules of algebra and 
equal powers of Q are collected, we obtain a polynomial in Q } ordinarily 
of degree n. The term independent of Q is simply the determinant | \a r8 \ |, 
which by hypothesis does not vanish. Now let F T8 denote the minor of 
Jrs in this determinant, taken with its proper sign. F rs is a polynomial 
in Q , ordinarily of degree n- 1. 
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Now operate on the first of (5) with F lm , on the second with F m , and 
so on, and add. Then 

fra ft* 88 F rm ( V r + QS r ) , (8) 

where summation with regard to s is understood on the left, and with 
regard to r on both sides. The operator acting on a particular x 8 is 
F rm f ra (r = 1, 2, ...ra). If s-m, this is equal to the determinant D\ 
if s 4= m, it is a determinant with two columns equal and therefore is zero. 
Hence (8) is equivalent to 

= Frm («r + QSr). (9) 

Now if all the S r are integrable within the range of values of t con- 
templated, the expression on the right of (9) is also integrable. Also since 
the function JD (z), obtained by replacing Q in JD by a number z, is 
regular and not zero at z = 0, the function l/D (z) is expressible as a 
power series in z with a finite radius of convergence. Define D ” 1 as the 
power series in Q obtained by putting $ for z in the series for l/D (z). 
Then operate on both sides of (9) with JD~\ We have 


D-'DXn, = D^Fm (Vr + QS r ) (10) 

But, since series of positive integral powers of Q can be multiplied ac- 
cording to the rules of algebra, D~ X D gives simply unity, and we have 
the solution 

X m = D- 1 Frm (tv + QSr) (1 1) 

This gives a complete formal solution of the problem. Its form is often 
convenient for actual computation, especially for small values of t; but 
it can also be expressed in finite terms. For the v r are constants, and 
(11) is therefore of the form 


_ fiui CQi) . F m 

J)(Q) + D(Q) 


QS r , 


( 12 ) 


where and F rm are polynomials, usually of degree one less than JD. 
Since the determinant | \a r(i \\ is not zero we may denote it by A , and D 
is the product of n linear factors, thus 


D~A (1 “ <*1 Q) (1 “ 0^(2) (1 — a n$)j (13) 

where the as will ordinarily be all different. Then <t>m(Q)l-D(Q) can be 
expressed as the sum of a number of partial fractions of the form 


L 

1 - o.Q 


(14) 


But by 1.12 (3) this is the same as Le at . The part of the solution arising 
from the u’s can therefore be expressed as a linear combination of expo- 
nentials. 
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The justification of this decomposition into partial fractions is that 
it is a purely algebraic process. Hence if the partial fractions and the 
original operator are all expanded in positive powers of Q, and like powers 
of Q are collected, the coefficients of a given power of Q will be the 
same in both expansions, and the equivalence is complete. 

Exceptional cases will occur if D contains no term in or if two 
or more of the « ? s are equal. If the term in Q n is absent the expansion 
of the operator in partial fractions will usually contain a constant term. 
If the term in Q^ 1 is also absent we must divide out, and the expansion 
in partial fractions will contain a term in Q ; this will give a term in t 
on interpretation. 

If several of the a’s are equal, the expression in partial fractions 
will involve terms of the form M(1 - aQ)~ r . These can be interpreted 
by direct expansion; but another method is more convenient. Starting 
with 

(l~a = ( 15 ) 

let us differentiate r - 1 times with regard to a. We find 


so that 


(r -1) ! Q r_1 _ fr-ljat 

(1-aQy 6 ’ 

(16) 

t T ~ J . 

(i-«Q) r (r-i)r 

(17) 


A rather different form of resolution into partial fractions from the 
ordinary one is therefore necessary if each fraction is to give a single 
term in the solution. Instead of having constants in all the numerators 
we must have powers of Q, the power of Q needed in any fraction being 
one less than the degree of the denominator. But if we write p~ l for Q 
the fraction in (16) is algebraically equivalent to p/(p - a) r . In this form 
the power of p in the numerator is independent of the degree of the 
denominator, and it is easier to resolve into partial fractions of this form 
than into those involving Q directly. 

The above remarks apply to the interpretation of the effect of the 
initial conditions ; that is, of the first term on the right of (11). To find 
the effect of the terms in the S ’ s we may expand the operators similarly. 
To interpret (1 - a^)- 1 QS in finite terms, we note that it is the solution 

of elf ~ ^ = ® ^at van ^ s * ies with *• This is easily found by the ordinary 
method to be 

y = e at Q(Se-*% 


which is the interpretation required. 


( 18 ) 
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To sum up, we can solve a family of equatio: 
first integrating each equation once from 0 to t 
ing for the initial conditions. This gives a set of eq 
(3). The subsequent process for deducing the operation 
is exactly the same as if the operators f r , were numbers, 





jd to tj allow- 


ancT the** 

ordinary rules of algebra were applied. The solution can be evaluated 
by expanding the operators in ascending powers of Q and evaluating 
term by term; this in general gives an infinite series. Alternatively it 
can be obtained by resolving into partial fractions and interpreting each 
fraction separately ; this gives an explicit solution in finite terms. 


1.4. Heaviside’s method is equivalent to that just given; it differs 
in using another notation, which is not quite so convenient in a formal 
proof of the theorem, but is rather more convenient for actual applica- 
tion. In Heaviside’s notation the operator above called Q is denoted by 
p~ l * At present we need not specify the meaning of positive powers of 
p ; negative integral powers are defined by induction, so that p~ n denotes 
Q n . We have noticed that the passage from 1.3 (5) to 1.3 (11) is a 
purely algebraic process. Consequently if all the equations 1.3 (5) were 
multiplied by constants before the algebraic solution the same answer 
would be obtained. Suppose then that we write the general equation 


1.3 (5) in the form 

(a rs + b r8 p ~ 1 ) x 8 = a r8 u a +p~ 1 S r (1) 

Multiply throughout by p as if this was a constant. We get 

( a r8 p + b ra ) x 8 = a r8 pn 8 + 8, (2) 


On solving the n equations of this form we shall obtain a solution 
identical with 1.3 (11) except that p~ l will appear for Q, and both 
numerator and denominator will be multiplied by the same power of p>. 
If the operators in the solution are expanded in negative powers of p , 
and p~ J is then interpreted as Q, the result will be identical with that 
already given. Comparing (2) with the original equations 1.3 (1) we see 
that the new form of our rule is as follows : 

Write p for djdt on the left of each equation; to the right of each 
equation add the result of dropping the b * s on the left and replacing 
the #’s by their initial values; solve the resulting equations (2) by 
algebra as if p was a number; and evaluate the result by expanding 
in negative powers of p and interpreting p~~ x as the operation of 
integrating from 0 to t. 

This is Heaviside’s rule. In what follows the equations (2) will usually 
be called the subsidiary equations. 
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To obtain the solution explicitly, we put e rs for a rs p + b ra , denote our 
determinant 


$11 

#12 

... 

$21 

$22 

• • • $2 n 

$m 

$«2 

• • • 


by A, and denote the minor of e r8 in this determinant, taken with its 
proper sign, by E rt . Then the solution is 

(***»«» + S r ) (3) 

Since the determinant formed by the a! s is not zero, A is of degree n in 
p, while i? m is at most of degree n - 1. The operators can therefore be 
expanded in negative powers of p> as we should expect; positive powers 
do not occur. All terms after the first vanish with t; the first is 

a rs u„ (4) 

where A ra is the minor of a rs in A . But 

ArmCtrs^A (m = S) | 

= 0 ( m*s)r W 

and (4) reduces to u mi as we should expect. This verifies that the solu- 
tion satisfies the initial conditions. 

1.5. To interpret the operational solution 1.4 (3) in finite terms we 
require rules for interpreting rational functions of p operating on unity 
and on other functions. We have already had the rules 

= Q> i?" 2 = Q * ; and soon, (1) 

p~ x l = p~ 2 l = = ~ ; and in general p~ n 1 = ~ . .,.(2) 

If unity is replaced by Heaviside's 'unit function,’ here denoted by 
H(t), which is zero for all negative values of t and 1 for all positive 
values, we still have 

p- n m=~ ( 3 ) 

when t is positive, but it will vanish when t is negative. We can also 
replace the lower limit of the integrations by — oo or by any negative 
constant without altering this interpretation. 
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Again, we shall have when t is positive 
_£_ = _J_ = e «t 

p-a 1-aQ ’ 

P _ O’ 1-1 _ tf" -1 at 

(p^T = (1 — aQ) n ~ (w-lj! fi “’ 


■(4) 

(5) 


1 = e“ e - 1, (6) 

jp-a jp-a 

where the function operated on may be either unity or H (t). In the 
latter case all the operators will give zero for negative values of t. 

The operators in 1.4 (3) are of the form f(p)/F(p ), where /Cp) and 
F{p) are polynomials in p, and f(p) is of the same or lower degree 
than F(p). If F(p) is of degree n it can be resolved into n linear 
factors of the form p — a. Then provided that the as are all different 
and none of them zero we have the algebraic identity 


whence 


Ap) ^ /(O) . ^ 

pF (p) pF (0) a «F' (a) p-a 9 

/GO -fM . 3 /(«) p 

F(pj~F(0) aF'(a)p-a- 


•(7) 


( 8 ) 


If this operates on unity or H (f) we have therefore for positive values 
of t 

/(?>), /(O) . X /(«) (Q) 

To justify this we notice as before that (8) is a. purely algebraic iden- 
tity, and therefore if both sides are expanded in negative powers of p , 
beginning with constant terms, the expansion of the two sides will be 
identical, and on interpretation in terms of integrations will give the 
same result. 

The formula (9) is usually known as Heaviside’s expansion theorem ; 
but as Heaviside’s methods involve two other expansion theorems* it 
will be called the ‘partial-fraction rule’ in the present work. 

If some of the a’s are equal or zero, the expression (7) considered as 
a function of p will have a multiple pole, and its expression in partial 
fractions will contain terms of the form (p - a)“ 8 , where s is an integer 
greater than unity, and a may be zero. Then f (jp)/F(jp) will contain 
terms of the forms or p/(p — a) 8 , which can be interpreted by 

means of (2) or (5). 


* Namely, expansion in powers of Q or and interpretation term by term; 
and expansion in powers of e~P\ where h is a constant, as in 4.2, 
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By means of these rules we can evaluate all the expressions for the 
part of a? m in 1.4 (3) that depends on the initial values of the x' s. If the 
S’s are constant for positive values of t, as they often are, the same 
rules will apply to the part of the solution depending on them. If they 
are exponential functions such as eP\ the easiest plan is usually to 
rewrite this as j pj(p - /a) and reinterpret. Thus 


Ap) 

F{p) 


. Ap) 


p 


( 10 ) 


AjA , t <i , a /(<* 

F(ri e +3 («- /i ) J P'(a) 


0 °* ( 11 ) 


If S is expressed as a linear combination of exponentials we can apply 
this rule to each separately. This is applicable to practically all func- 
tions known to physics. 

Alternatively we can resolve the operator acting on 8 into partial 
fractions and interpret % r n 8 by integration and other fractions by the 
rule 1.3 (18) 

— fse-^dt. ( 12 ) 

P-clJq v 


This completes our rules for solving a set of linear equations of the 
first order with constant coefficients. In comparison with the ordinary 
method, we notice that the rules are direct and lead immediately to a 
solution involving operators, which can then be evaluated completely 
by known rules. If it happens that we only require the variation of one 
unknown explicitly, we need not interpret the solutions for the others. 
In the ordinary method we have to find a complementary function and 
a particular integral separately. To find the former we assume a solu- 
tion of the form x 8 -X s e ai and on substituting in the differential equa- 
tions with the 8’s omitted we find an equation of consistency to 
determine the n possible values of a, and the ratios of the Vs corre- 
sponding to each. The particular integral is then found by some method, 
but it does not as a rule vanish with t. The actual values of the Vs 
are still undetermined, and the value associated with each a must be 
found by substituting in the initial conditions and again solving a set 
of n simultaneous equations. The labour of finding the equation for the 
a’s and the ratios of the Vs corresponding to one of them is about the 
same as that of finding the operational solution of Heaviside's method; 
the rest of the work is avoided by the operational method. Further, if 
some of the Vs are equal or zero considerable complications are intro- 
duced into the ordinary method, but not into the operational one. 
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1.6. The above work is applicable to all cases where A is not zero. 
If it is zero, we can show that there is some defect in the specification 
of the system. A system is adequately specified if when we know the 
values of the dependent variables and the external disturbances at any 
instant the rates of change of the dependent variables are all determinate. 
Now if A is zero, let us multiply the rth equation 

+&,)*. = & ( 1 ) 

by A rm , the minor of a rm in A , and add. Then 

4 ™ a n ^ = Am, (, S r - b rs x,) ( 2 ) 

The coefficient of dxjdt on the left is A rm a r8 . If s^=m, this is a deter- 
minant with two columns equal, and is therefore zero ; if s = m, it is equal 
to A, and is again zero. Thus the left side vanishes identically, and the 
x's are permanently connected by the relation 

A rm (b r8 x 8 -Sr) = 0 (3) 

If the coefficients A rtn b r8 vanish, there must be a permanent relation 
between the JS r , and one of the original equations is a mere logical 
consequence of the others ; then we have not enough equations to deter- 
mine the derivatives. If they do not vanish, the values of the x s when 
t~ 0 are connected by a definite relation, so that the initial values of 
the x s cannot be assigned independently. 

1.7. The method is most easily extended to equations of higher 
order by breaking them up into equations of the first order. Thus if we 
have an equation of the second order such as 

d*x dx . , « 

dt +a dt +b ‘' B ~ S > CO 

we introduce a new variable y given by 

dx .. \ 

dt~y =0 ’ (2) 

and the original equation can be replaced by 

^-+ay + bx=S. (3) 

We have now two equations of the first order in x and y. If initially 
x = x 0 and dxjdt = , the subsidiary equations are 

px-y=px 0) (4) 

(p + a)y+bx=px i +S. (5) 
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Solving by algebra we find 

_ ( p a + ap ) # 0 +pa i + 8 
p 2 + ap + b 

To interpret, put 

p 2 + ap + b = (p - a) (p — /?), 
and apply the partial-fraction rule. We find 

# = -~p fa- fo) ^ Oi - cur«) *** 


.( 6 ) 

‘(7) 


~ [«■* J* Se-'t at - &' J* 8r+* *], ...( 8 ) 


which is easily shown to satisfy all the conditions and therefore to be 
the solution required. 

1.71. A few illustrative examples may be given, ( p is written for 
djdt from the start.) 

1. (p 3 + 4p + 3)# = 1; = # x = -2. 

We consider the subsidiary equation 

( p* + 4p + 3 ) # = 1 — 2p + 3 (p 2 + 4p) 

~ Sp 2 + lOp 4 - 1 . 

_ ffi+lQp + l 
(p + l)(p + 3) 

= & + 30- t -£e-” 

on interpreting by the partial-fraction rule. 

2. (jp 2 +5jt? + 6)#=12; # 0 = 2; #i = 0. 

Consider (p 2 + 5p + 6) # = 12 + 2 ( p* + 5p ). 

Cancelling the common factor, 


Then 


# = ; 


3. 

This can be written 


#= 2 . 

(p + 3)« = e- a< ; y„ = 0. 

(y + S^-^i 


or 


#= ; 


p 


4 . 


Here 




0 + 2)0 + 3) 
(j» a + w a ) ir=0. 

_^£o+££i. 


£r + 5 


= #0 COS 


sm 
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5. (^ + 2) a # = tf a e~ 2 * ; ^ O =0; ^ = 0. 

The expression on the right is equivalent to the operator %p/(p + 2) 3 . 


Hence 


2 p 


4! 


~ — t* a- * 
12* 9 • 


{p^r 

We notice the advantage of Heaviside’s method in avoiding the use 
of simultaneous equations to determine the so-called ‘arbitrary ’constants 
of the ordinary method. In particular in the second example the data 
have a property leading to a simple solution. Heaviside’s method seizes 
upon this immediately and gives the solution in one line; with the 
ordinary method the simultaneous equations for the constants would 
have to be solved as usual. 

In general we may say that the more specific the problem the greater 
will be the convenience of the operational method. 


1.72. The method just given can be extended easily to a family of 
equations of the second order. If the typical equation is again written 

&rs X a ” $r j (l) 

where now e ra = a rs ^ + b ra j t + c ra , (2) 

we can introduce n new variables y u y a , y n given by 

3?-*-° <«> 

thus treating the first derivatives of the #’s as a set of new variables. 
Then the equations (1) are equivalent to 

Or. + ( b ra y , + c ra x a ) = 8 r , (4.) 

and (3) and (4) constitute a set of 2 n equations of the first order. Suppose 
also that when t = 0, ® a = u a , ys-%* According to our rule of 1,5 we 
must replace djdt by p and add pu a to the right of (3), and a ra pv a 
to the right of (4). We have now to solve by algebra, and may begin by 


writing the revised form of (3) 

y 3 =P (x 6 -u a ) (5) 

When we substitute in the modified form of (4) we get 

a ra p 2 (cc a - u a ) + b ra p (iv a - u 6 ) + c ra x 8 = Sr + a ra pv a , (6) 

or, on rearranging, 

(a ra p* + b ra p + c ra ) X a = (a r8 p* + b r$ p) u 8 + a ra pv 8 + &, (7) 


m 
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We have thus n equations for the #'s to solve by algebra, and the 
solution can be interpreted by the rules. 

In Bromwich's paper ‘Normal Coordinates in Dynamical Systems' a 
method equivalent to the operational one is applied directly to a set of 
second order equations. First order equations, however, arise in some 
problems of physical interest, and merit a direct discussion. This, as we 
have seen, is easily generalized to equations of the second order. 

1.8. In all the problems considered so far the operators that occur in 
the solutions are expansible in positive powers of Q, or in negative powers 
of p, and we have seen that so long as the series obtained by replacing 
Q by a number have a finite radius of convergence, however small, the 
operational solution is intelligible in terms of the definitions we have 
had. But we have not defined p as such, because we have only needed 
to define its negative powers, and p cannot be expanded in terms of its 
own negative powers. Then has p any meaning of its own? Since it 
replaces djdt when we form the subsidiary equation we may naturally 
suppose that p means djdt , and this is the meaning sometimes actually 
attributed to it; but care is needed. We recall that when the subsidiary 
equation is formed a term like px Q appears on the right; but dx^jdt is 
zero, so that if we pushed this interpretation too far we should be faced 
with the alarming result that the solutions of the equations do not 
depend on their initial values. The fact is that though the operators 
djdt and Q both satisfy the laws of algebra and are freely commutative 
with constants, they are not as a rule commutative with each other. 
Thus 

!«/«-! I'/m-m, <» 

« % m - jT r « & -/(<> -/(o) <2) 

Thus the operators p and Q are commutative if, and only if, the function 
operated on vanishes with t. It appears from (1) that djdt undoes the 
operation Q if djdt acts after Q. With this convention we can identify 
p , the inverse of Q, with djdt When p and § both occur in an operator, 
the Q operations must be carried out before the differentiations*. 

My own view is that the method of expansion in positive powers of 
D, which is given in text-books of differential equations for finding 
particular integrals, should be abandoned because it leads to intelligible 
results only in two very special cases, and should be replaced by the 

* Of. Heaviside, Electromagnetic Theory , 2, 298, 
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preseat one. The treatment of such equations by the present method 
is as follows. Suppose that the differential equation is 
■d\ 

I /)! ^ P 

dt n 


♦(*) 




d n ~ l y , 

+ +- + a»3f-f<fa 


where f(t) is known, and that when t- 0 

dy 

y=y»\ T+ = yi> 


d n ~ l y_ 

— y»-i • 


We integrate n times from 0 to ti 3 determining the constant of integration 
at each stage. The resulting equation is 

1 + &i Qy +...+% Q n y = (1 + ct x Q + ... + c^n-i Q n ~ 1 ) 

+ (Q + «i Q 2 + . . . + ^_ 2 Q 71 - 1 ) y l 

+ 

+&-'%-!+ or /(t). 

We write Q as p~ l and write the subsidiary equation as 

<t>(p)y = (p n + a iP n ~ 1 + — +a„)g=(j> n + a 1 p n - 1 + ... + o a - 1 p)g 0 

+ (p n ~ 1 + a 1 p n - a + . . . + ®»-ap) y t 


+P!/n-l +/(*) 

. , , . . +/00 

say. Then the operational solution is 

The first part expresses the effect of the initial conditions, and in the 
usual nomenclature gives the complementary function. The part depend- 
ing on /(t) gives a particular integral, which can always be obtained 
formally by 1.5 (12). 

The text-book method (due to Boole) is applicable if f(t) is a poly- 
nomial, or if it is expressible in the form e^F (tf), where F(i) is a 
polynomial. Suppose first that our solution is 


where m is an integer. We write this as 
ml 


ml 


y 4>(p)p m (p - <0(p- **») ---(p - Op m ‘ 

The operator must be broken up into partial fractions. If all the a’s are 
distinct and different from zero, the fractions arising from them are 
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-which would in Boole’s treatment be part of the complementary function. 
The contribution from the part near p = 0 is to be found by expanding 

V -»*• _J_ 

* P m '*{p) 

in ascending powers of p until we reach a constant term. Then 

Y _ ml 1 \ 
p m a n + ctn^ip + a n \ / 

ml /t m r- 1 f \ 

a n \m\ a n (m- 1)1 "7 

a<n Cln 

taken as far as the .constant term. But this is precisely the result given 
by Boole’s method when we expand 1 /<j> (D) in ascending powers of D 
as far as D m and make it operate on t m . To this extent therefore the 
methods are equivalent; but Heaviside’s method gives the integral that 
vanishes with its first n - 1 derivatives at t = 0, while Boole’s leaves the 
coefficients of the exponentials to be determined separately. 

Now consider the function p~'[e!* F(t)\ where F(t) is a known in- 
tegrable function. We have 


p 1 F(t)J = jf F(t) dt = jf <Fd {p~* F (t)} 

= fetp-'Fit) J‘ - pJ‘<Fp~'F(t) dt 

= tf* [p-i _ ap ~* + a * p - » _ . ]F(tl 

aU integrals vanishing at # = 0 since F(t) is integrable. Thus 
p- 1 [e* F(t)] = e« f ~L Fit'), 

•n . P + a v 7 

Kepeating the process, we have 


p-*[0*F(t)] = e ai - 


.^(0, 


(p + a) 7H ' 

“ d lWor *' if V* M >■« expansible ip negative power, otp, 

Tb Si^lZT'“ i “ g “f ’ " k «“ oeStaMy tortaont. 

given in textbSL of daS f *■«* tan tlioso 

fferential equations for the dctormination of 
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particular integrals. For instance, we have interpreted l/(p - a) as 
- (e at -l). In the ordinary method, we expand this operator in ascending 
powers of p and then interpret p as d/dt Thus 


+ i (3) 

p—a \a or / a 

This assumes that the sum of all the terms on the right after the first 
is zero. But 

( — +-)i= -/-' I w 

\p-a a J 

By our rules we should interpret the function on the right as ~ e at . If 


we write itas-.jo. — — . 1 and interpret the first p as d/dt, it gives 

a * p - a 


a p — a a dt [a ' ') a 


1 1 


■(5) 


agreeing with the former interpretation. But if we write it as - 

and interpret the p as meaning d/dt, it is clearly zei’o. The difference 
arises from the fact that the operators d/dt and (p - a)” 1 are noil- 
commutative ; the right result is the one obtained by operating with 
(p - a)" 1 first 

We may notice, incidentally, that whereas the series in powers of Q 
always give convergent series on interpretation, the same is not true 
of the corresponding series in powers of d/dt For instance, if S = t r f 

where r is fractional, the series for . S in ascending powers of p 

P 

diverges like Apart from the greater internal consistency of 

Heaviside's methods, then, they are capable of much wider application. 


1.9- An interpretation of the operator e hp is suggested by Heaviside, 
as follows : 

•*/(*) = (l+Ap + §f + -)/00 

-/(*) + V (*) + £/" (*) + ••• 

=f(t + h) ‘ (6) 

by Taylor's theorem, if we interpret^ as d/dt The operator then increases 
the argument of the function by h. But this is clearly dangerous ; for 
suppose that / (jf) is 0 for - 1 < t < 1, and equal to t* — 1 for t outside 
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these limits. If in the series on the right we put t = 0, every term 
vanishes, and apparently for all values of h 

f Qi) =/( 0 ) + hf ( 0 ) + ... 

= 0 00 

But/ (A) is not zero if | h | exceeds 1. The fact is, of course, that Taylor’s 
theorem is not applicable to a function of this type, but only to functions 
expressible by a single power series throughout their range of existence. 
It actually appears, however, as we shall see later, that the formula (6) 
is more general than Taylor’s theorem and that its justification does 
not depend on the interpretation of p as djdt . At present we may simply 
take it as the definition of the operator; we define e hp f (t) to mean 
f(t + h\ We must now examine whether this operator commutes with 
Q and with djdt . We have 

Qf(t) = * ' j*f(t ) dt = dt, (8) 

rt rt+h 

Q e hv f(t) = Q/(t + A) = Jy(t + /i)dt= j h f(t)dt, ...(9) 

rh 

so that the operators d* and Q are commutative if j f(t) dt = 0. This 

restriction is less serious than might appear, because many of our 
functions vanish for all negative values of t. Then if h is negative the 


condition is satisfied. 

Also + ..(10) 

it = iA* + *) CD 

so that e hp and djdt are always commutative. 


The operator e hp is useful in enabling us to find a theorem that plays 
in Heaviside’s method a part closely analogous to that played by 
Fourier’s integral theorem in the ordinary method. Instead of treating 
trigonometric functions as fundamental, Pleaviside uses his 'unit 
function,' here called H(t), which is unity for positive and zero for 
negative values of t. Then 

e hp ff(t) - H(t + h) (12) 

and is equal to 1 if t > - h and to 0 if t < - k Then if h x < 

(er*h-e-p\) H(fi) = 

= 1 if h 1 <t< h 
= 0 if t < k 
and = 0 if t > Ji 2 



( 13 ) 
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If, further, h<h i <h ! ><...< h n , 

{(e~*h - f(1 h ) + (fr&% - e~*h) f(/i 2 ) + . . . 

+ - e-* h *)f(]in- 0 } #(*) 

= 0 for £ < /^/(Ai) for A x < and so on. ...(14) 

If then the subdivisions of the interval hi to li n become indefinitely 
numerous and close, we can approximate as closely as we like to a 
function f (t) by a sum (e~ ph r - e~ ph r+ 1 ) In the limit we 

can replace this sum by an integral and write 

/(<) = - [ = _J(h)d[e~*H(t)] ( 15 ) 

Ilae function under the d sign being discontinuous at h = t, the integral 
is of Stielljes type*, and its meaning is really defined as the limit of (14). 
But we can often formally integrate with regard to h as if p was a 
number, and obtain a function of p operating on II (t). In this way a 
given function of t can be expressed in the operational form, thus 

m=g{p)H(t) (16) 

This device often simplifies the evaluation of the result of an opera- 
tion on a function ; for if our solution is 4> (p~) f if), say, we can interpret 
it as <j> (p) g (p) II ( t ), and then interpret directly the result of operating 
on //(if) with <t> (p) g (p). 

If in the equations 1.3 (1) all the S’ s are zero except one, and that 
one is //(if), it gives a contribution to any coordinate of the form 

n:=<i>(p)H(t) (17) 

Suppose now that we replace H (t) by H{t-h). The corresponding 
contribution to the solution is <j> (p) II (t- A). We can now extend (17) 


to a general S, for if, as in (14), S=f(h,) from t = h r to h r+1 , we can 
write the contribution to the solution . 

x (p)f(h) {II \t -hr) - II(t - Ar+0}, (18) 

and by superposing contributions from different intervals of time and 
proceeding to the limit 

x = ~ f <l>(j))S(/i) dll (; t-h ), ( 19 ) 


which is in general a Stieltjes integral. In this way we obtain a solu- 
tion, with II (t) regarded as the fundamental function, which replaces 
decomposition into trigonometrical or exponential functions as in the 
ordinary treatment. 


Hobson, TJicory of Functions of a Meal Variable , 1 (1921), 507. 
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CHAPTER II 

COMPLEX THEORY 

2.1. All the interpretations found so far for the results of operations 
on unity or on H (t) are special cases of two closely related general rules, 


as follows. If $ is an analytic function, 

m 

(a) 


In the former definition the integration extends around a large circle in 
the complex A plane. In the latter the integration is along a curve from 
c - too to c + too , where c is positive and finite, and the path is such that 
all singularities of the integrand are on the left side of the path (that is, 
the side including - qo). These rules are both due to Bromwich*. We 
have to show that they agree with the interpretations already given. 

Considering first (1) we have, when n is a positive or zero integer, 



and the only pole of the integrand is at the origin, where the residue is 
t n /nl Hence when n is a positive integer or zero 


agreeing with the rule of the last chapter. 

If <£ ( p ) is expressible as a power series in such that on replacing 
p- 1 by a number we get a series with a finite radius of convergence, we 


can write 

<f> (p) = a Q + a x p^ + ... + a n p~ n + ..., (5) 

and by the rule of the last chapter 

t n 

$ (p) . 1 =Oo + <ht + • •• + «» + (6) 

7l I 


But by (1) the interpretation is 

« 

* ‘Normal Coordinates in Dynamical Systems/ Proc. Lond . Math. Soc 15 (1916), 
401-448. 
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which, we see from (3) and (4), agrees with (6) provided term by term 
integration is justifiable. To justify it, we notice that the series in 1/A. 
converges to a finite sum whenever |l/A| is less than a certain constant 
r; hence if we write 

(A) - a Q + ~ + Rn (A) (8) 

and |A| > 1 fr, we can choose n large enough to make 

\%nW\<* ( 9 ) 

everywhere on C, however small € may be. Hence 

$ (p ) . 1 = a 0 + a x t + ... + <v- 2 ^ + cfrj c W ^ (10) 

If then G is a circle of radius 1 jr or more, the last term is less in ab- 
solute value than e tlr e. Hence if r is chosen to begin with, we can take 
e small enough to make e t/r e less than an arbitrary small quantity rj, 
and then n large enough to make | R n (A) | less than this e. Thus 
(p) . 1 differs from the sum of the first n terms of the series by 
less than an arbitrarily small quantity rj, and therefore term by term 
integration is legitimate. 


2.11. Some of our interpretations in finite terms can easily be proved 
directly. Thus 

p . __** ,, 1 at f ^ 

(p-a) 71 ' (A — a) 11 2m 6 Jq( k — a) n } 

which has a pole of residue jf n_1 /(ra - 1) ! at A = a. Hence 

t n-l 


P 


1 = 


V, 6 at 


(p-af (fl — 1)! 

as in 1.5(5). 

'Also we can verify the partial-fraction rule. Thus 


./ (p) , = JL (H 
F(p ) ' 2tti J a 


fW 


A(\-aj (A — Oj,) ... ( A- <%) 


dX 


\ T ** - Ut f / 


F( 0) ' ~aF'(a)' 
on evaluating the residues at the poles 0, « 1; a a , ... a,,. 


•CD 


■(2) 


2.12. We can also prove directly that 1.4 (3), when interpreted by the 
rule (1), is the correct solution of equations 1.3 (1). If we denote by 
e rs (A), A (A), and E rs (A.) the results of replacing p by A in e r8i A, and 
F r8) the part of 1.4 (3) not involving the S r is equivalent to 




1 [ JBm(K 

2m] o A (A) 


v r e Kt dK 


a) 
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instituting in the differential equations and differentiating with regard 
* under the integral sign, we get for the left side of the £th equation 

euX * = e,a M Vr eKtdX " (2) 


»(A )v r & t d\, (2) 


lt -EnWeic (A) = 0 if r * Tc } 

, , =A(A)if r = k] W 

d therefore 

= -(*) 

lus the differential equations are sa tisfie d 

Now considering the initial conditions, we put t = 0, and find 

■« 

t v r = a ra u 8 = ^ {e ra (A.) — b ra ] u 8 ... (6) 

1 Xm= 2 mj c A™ (A.) M ^ “* 

= 2^J 0 Wm ~AA(A) ^ ( 7 ) 

w A (A) is a polynomial of the rath degree in A, and E m (A) one of 
i (»-l)th degree. When |A| is great enough the second term in 
integrand is therefore of order A -2 , and gives zero on integration. 


I the initial conditions are satisfied. 

tfext, suppose that instead of the S’ a being zero they are exponential 
ns of the form iV‘. Since the solutions are additive it is enough 
suppose that the initial values of the y’a are all zero, and to consider 

effect of an exponential term in the first equation alone. Then the 
ations are 

e nVi + <hai/a + . . . + e ln y K = Pe* 1 = -Si 

„ , N P-F- ■ (9) 

0, •«?/„ = 0 (r 4= 1) ' 

i operational solution is 


JUn Pp 
A p-ij.' 


is to be interpreted as 

„_J ( Bu{ A) P 

“ 2m J 0 A (A) A — p 6 ^ (If) 


4 '* 
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Substituting in the differential equations we find that the left sides 
of all vanish except the first, in consequence of the integrand con- 
taining as a factor a determinant with two rows identical; the first 
gives 

~ f ~dk=Pe^. ( 12 ) 

2m JgA'-P' 

The differential equations are therefore satisfied. 

When £ = 0 and |A| is large, the integrands in (11) are all 0 (A.” 2 ) and 
the integrals are therefore zero. Thus all the #’s vanish with t 

Since all functions that occur in physics can be expressed by Fourier’s 
theorems as linear combinations of exponentials, the proof that the 
solution can be carried out by operational methods is complete. We 
shall see later, however, that Fourier’s theorems are not so convenient 
to use as the formula 1.9 (19) based on the function H it). 

2 . 2 . Now let us consider the integral 2. 1 (2), taken along the path L . 
If we suppose the contour completed by a large semicircle from c + too 
to c- ico by way of -co , all the singularities of the integrand are within 
the contour, and therefore the integral around it is the same as the 
integral (1) around a large circle. Now if t is positive, and (A)/\ tends 
to zero uniformly with regard to arg A as A tends to co , the integral 
around a large semicircle tends to zero as the radius becomes indefinitely 
large, by Jordan’s lemma *. Hence the integral along L is equivalent to 
the integral around a large circle. Thus if t is positive and 

$ (A)/ A = 0(A -W ) 

when A is great and n is positive, (j;) II (t) is the same as <j> ( p ) 1. 
But if t is negative, the integral around the large semicircle on the 
negative side of the imaginary axis is no longer an instance of Jordan’s 
lemma, and this result no longer holds. The integral around a large 
semicircle on the positive side of the imaginary axis is, however, re- 
ducible to the integral along L, since by construction there is no 
singularity between these two paths. If then <j> (A)/A = 0 (A“") when A is 
great, the limit of the integral around this large semicircle is zero by 
Jordan’s lemma, aud therefore </> (p) II (t) is zero. 

Thus if (p) is expansible in descending powers of p, beginning with 
a constant or a negative power of p, we shall have 

<£ ( p ) II{t) ~<}>(p)l when t > 0 \ 

<f> {p) H{t) = 0 when tf < 0 J ’ ^ 

and the results agree with those of Chapter I. 

* Wilittaker and Watson, Modem Analysis, 1915, 115. 
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2.3. The form 2.1 (1) suggests an interpretation of p n 3 where n is 


a positive integer. For it gives ) 

p n .l=^-f X'-'eX'dX 

* ViriJ o 

= o, (1) 


since the integrand has no singularity in the finite part of the plane. 
This result is the same as is obtained by interpreting p as d/dt. 

If <f> (p) is expansible in negative powers of p as in 2.1 (5), we have 

S<* «■!)-: 

-h,L^ mdK - ® 

since differentiation of a contour integral may be carried out under the 
integral sign. 

But (2) is the interpretation, according to 2.1 (1), of p <j> (jp ) . 1. Hence 
if we have an operator involving both definite integrations and differen- 
tiations, and we replace definite integration by p~ x and differentiation 
by p, and interpret according to 2.1 (1), we get the result of carrying 
out the differentiations after the integrations. The complex integral 
interpretation is therefore consistent with 1.8. 

2.31. Again, if the operator <£ (p) acting on unity gives /(£), and we 
seek for an interpretation of <? v f (0> we have 

and the residues given by the various terms are the same as in the 
interpretation of $ (p) . 1 except that t is replaced by t + h. Hence 

( 2 ) 

If the operator ( p ) acting on H(t) givesf(t), we may proceed similarly, 
using the line L instead of C; again the form (2) holds. This gives the 
justification for 1,9 (6),- for any bounded function can be expressed in 
terms of H{i) by 1.9 (15), whether it is expansible in a power series or 
not. We had 

/(*)=- P (8) 

J ft, =1—00 

and the right side, interpreted by 2.1 (2), gives 

-GiU/n*”**- w 
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On L we may put 

k = c + IK, 

where c is real and positive. Then 


.(5) 


= ~ J°° f " fQije^e^-^dhdK 

= ~J j f(h) +* cos K {t—h) dhdK. (6) 


We may rewrite this as 

/(*) r« = - f " f“‘/(A) e- rf ‘ cos k (# - A) (7) 

TT J -co JO 

so that (G) is equivalent to Fourier’s integral theorem applied to the 
function /(#) <T C *. It is therefore more general than this theorem applied 
to f(t) directly j for with c equal to zero the integrals may diverge, 
but the theorem may be saved by giving c a finite positive value. But 
just as the interpretation of is usually independent of the 

value of c, so here the value of the integral is independent of c so long 
as it exceeds a certain value, which in many practical cases is zero. 

In practice the form 2.1 (2) or a modification of it is generally used 
in preference to 2.1 (1). In problems involving a finite number of 
ordinary differential equations of the first order the forms are equivalent 
for all positive values of the independent variable ; and since the inde- 
pendent variable is usually the time, and we require to know how 
a system will behave if it is in a known state at time 0 and is afterwards 
acted upon by known disturbances, it is as a rule only positive values 
of the time that concern us. Then an integral along L is equivalent to 
one around G. Indeed if instead of L we take a curved path beginning 
in the third quadrant, passing the origin on the positive side, and 
receding to infinity anywhere in the second quadrant, provided we still 
keep all poles on the left of the path, we shall obtain the same value 
for positive values of t ; this integral in general will diverge for negative 
values of t. 

But when we come to deal with continuous systems it usually happens 
that when the operational solution is interpreted as an integral, 
the integrand has an infinite number of poles, and that no circle however 
large can include all of them. Hence the contour integral 2.1 (1) cannot 
be formed. But the line integral along L usually still exists. Again, it 
may happen that <j>(k) has a branch-point at the origin or on the 
negative side of the imaginary axis. Here again the contour integral 
does not exist, but the line integral does. 
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Consequently in most of the writings of Heaviside and Bromwich, 
when a function of an operator <£ ( jt?) occurs without the operand being 
specified explicitly, it is to be understood that the operand and 

that the interpretation 2.1 (2) or a suitable modification of it is to be 
adopted. This rule will be followed in this work when we come to deal 
with continuous systems ; and if it is also supposed to be adopted in the 
problems of the next chapter no harm will be done. 

Considerable latitude is seen to be permissible in choosing the path L, 
Since by construction the integrand is regular at all points to the right 
of L , we can replace L by any path on its positive side provided it ends 
at c ±ioo. Again, the systems we discuss are usually stable systems: 
that is, the poles of </> (A) are all on the imaginary axis or to the left of 
it. Then L can be taken as near as we like to the imaginary axis. L must 
not cross the imaginary axis if there are poles on the latter with 



Mg. 1. 1 

indefinitely great values of the argument. Bromwich's usual device is 
to take L parallel to the imaginary axis and any constant distance o 
from it; by our restriction that all singularities must be to the left of X, 
however, we avoid a possible ambiguity for unstable systems. By allowing 
L to begin and end within the third and second quadrants we ensure 
convergence in some cases when it would otherwise fail, but this is 
never, so far as I know, necessary in the actual solution of a physical 
problem, though it is often useful in the course of evaluation, as in 
2 . 4 ( 22 ).' 
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2.4. The principal operators involving branch-points are fractional 
powers of^. We require then an interpretation of p™ 11(f), where m is 
a positive or negative fraction*. By our rule 

(i) 

when t is positive. The integral converges if m < 1; but if L is replaced 
by a path L such that the ends are in the second and third quadrants 
and make finite angles with the imaginary axis, the integral will converge 
for any value of m. There is a branch-point at A. = 0. 

A contour containing L as part of itself, and such that the integral 
is regular within it, is as shown. Evidently the large quadrants make 
no contribution, subject to the proviso that if m > 1 we replace L by L\ 
The integral along L or 1! is therefore in the limit equal and opposite 
to one along CDEF \ If further m is positive the contribution from DE 
tends to zero with the radius of DE. Now on CD and EF 


\ = fxe Llr and pe~ tn , (2) 

respectively, where p is real and positive. Hence on CD 

XP-'ePdk = - = ^e^e^d^ (3) 

and on EF 


X m ~ l e xt d\ ( 4 ) 

Thus 

r° r 00 

2mp m n(f) = - J fjL in ~ 1 e" fit e mL7r dfJi — J f n ~ 1 6~ flt e~ miir dfA 

= 2i sin mir f /a” 1 ” 1 er^ dft (5) 

Jo 

But if we introduce Gauss’s II-function, defined for positive values of 
m by the equation 

II ( m ) = r (m -i- 1) = [ e^ l u m du 3 (6) 

Jo 

and for negative values of m by the inductive equation 

H (m) no-i) (7) 

we have J g = (8) 

Hence p m H(t ) = — ( w ) s i n mv ($) 

mir 

* Of. Bromwich, Proc. Camb. Phil. Sob. , 20 (19 

siilwi Llb ' B:|or; 
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Bat by a known theorem 

n (m) n (- m) = wztt cosec m?r, (10) 

tr m 


and therefore 




.( 11 ) 


n(-m) 

This has been proved for positive values of m. But if we define p ml II (t) 
for any value of m by 2.1 (2) with LI for L we have , 


p m+l = _L j^^t dX 

-LLk km *<n 


■ [t *”] ' aSJ w 

Tlie integrated part' vanishes at both limits. The integral is - - p m 11(f). 
When m is positive we can apply (11) and get 

t~m t - m -i 


p m+1 H(t) = 


m 


.(13) 


n(- w -i)’ 

which could also he obtained by putting m + 1 for m in (11). But wo can 
also put 

p m S(t)=- — p m+ 1 E(t) t ~ m ~ 1 t ~ m c,n 

which extends (11) to values of m that make m + 1 positive. Hence (11) 
js true for *i>-l. By repetition we can show that it holds for all 
values of m. If we change m into -n we have for all values of n, 

(i#) 

Since when «t is a positive integer U( n )=n\ and when n is a negative 
integer n (n) is infinite, the interpretations already adopted for integral 
powers of p are special cases of (15). “ 

If we differentiate (15) with regard to n, we have formally 

., . n (w) in (jj)p * 

and if then n tends to 0, W 1 '' 

- iogi?=io g «-n'(o)/n(o). 

But - H'(0)/n(0) is Euler’s constant y, and 
log p =mm y log t. 
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We notice also that 

_£!‘ _ = _ m tr-n-1 _ t -™- 1 _ _„ + i rr (A 

dt* •° U rfif n (— »a) n(-»ra) "w* 

( 16 ) 

«) - ajrsQ-y-JCO. ...(17) 

provided m is negative. Hence the differentiation is equivalent to 
multiplying by p and reinterpreting ; and integration is equivalent to 
dividing by p and reinterpreting, provided the result of integration 
vanishes when t = 0. 

In particular, since 

n(-i) a .r(j) = v, (is) 

we have 

pia(,> jkr ?'-* * "(«) 

and so on. Also 

PO) 

A related function that occurs in problems of heat conduction is 

where a is a constant, and q denotes p^. By 2.1 (2) 

r-BW-kL — *- < 8 » 

On L the argument of A- is between ± ^7r. Hence if a is positive the 
integral is convergent. 

Immediate expansion of <3~ oA * in a power series and integration term 
by term would not be legitimate, because all the resulting integrals 
after the first two would diverge when taken along the path L. But if 
instead we use a path L\ (21) is unaffected, and on L we can proceed 
in this way. Thus 

k it -if +•■•]*■ •••<»> 

All the positive integral powers of A give zero on integration. The other 
terms are equivalent to 

i-ap 

* 3! 5! 

2 r « i / « y . 1 ( a v i_i r a 

s/irl2 Jt S\2jt) 2 1 5 \2 Jt) •"/ 1 ert 2 Jt’ 

(23) 
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9 f w 

erf w — j q (24) 

By differentiation with regard to a we find 

^' aaff ^ = ^) e ~ mt - < 25 > 


We shall sometimes need an asymptotic approximation to erf w when 
w is great. We have 

l-erfw=-7“ f e~*dt=-i- f e~ u u~*du 

sjTT J w J wfl 

1 -wa -iTt 1 -a 1-3 -4 1.3.5 . -I 

= W \}~Z W + 2T2 W 27272 W + -J 

(-1)“ 1.3.5... (2 n - 1) T , , . 

+ V^ 2.2.2... 2 - JJ ( 2G ) 


on successive integrations by parts. But the last integral 

= 2 f e- tl t~ m dt<2e~ wi [ t~ in dt 
J W J W 


< W * w - 1 >, 

2ra- 1 ' 


•(27) 


so that its contribution to the function is less than the previous term. 
We have therefore the asymptotic approximation 


(x — erf vi) ~ -4- ["" w~ l — 4 w~ s + 

<sAL 2 2.2 


w~ 


1.3.5 

2 . 2.2 


w~ 




2.5. Carson’s method. In this method p is not regarded as an 
operator but as a quantity. If we return to 1.9 (15) 

/(#) = “ /*__'/(*) d [«- vh H{tj\ (1) 

we can write formally 

= f_"Pf(h)e-» h dh ( 2 ) 

and then 


g(j>)H(t)=f(t) (8) 

• /-/n q ^ ^ 1 1S k nown J 18 an integral equation to deter- 

r-i®, 8 re “5 nct,on “ ■'»’> " i»pBed by the oom.~o of 

gr , ra practice f(h) vanishes when h is negative, so that li is 


COMPLEX THEORY 33 

only allowed to proceed from 0 to go , and p is taken to have a positive 
real part. Then (2) becomes 

g (p) = J o e-* h dh, (4) 

an equation to determine /(h) for positive values of k. This is Carson’s 
integral equation. Here, however, (2) will be used. 

Now the integral 

, i fj c m mjp 

=/(0 (5) 

for the double integral in (5) is simply that in 2.31 (4) with p written 
for A.. Thus the solution of (2) is 

» 

This proves that Bromwich’s integral is the solution of Carson’s equation, 
as has been pointed out by Ldvy and March. 


CHAPTER III 

PHYSICAL APPLICATIONS: ONE INDEPENDENT VARIABLE 

3.1. An electric circuit contains a cell, a condenser, and a coil with 
self-induction and resistance. Initially the circuit is open. It is suddenly 
completed ; find how the charge on the plates varies with the time. 

Let y be the charge on the plates, t the time, K the capacity of the 
condenser, L the self-induction, and B the resistance of the circuit; 
and let E be the electromotive force of the cell. Write p for djdt. 

The current in the circuit is y, and the charging of the plates of the 
condenser produces a potential difference yjK tending to oppose the 
original e.m.f. Then y satisfies the differential equation 

j T-£=lj/ + JSy (1) 

Initially y and y, the current, are zero. Hence the subsidiary equation 
is simply 

(ip^ + Itp + j?) y = E, 


( 2 ) 
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and the solution is 


y= 


E 


LfP + Rp+jj? 


( 3 ) 


If now the denominator is expressed in the form L(p + a)(p + /3), the 
interpretation is, by 1.5 (9), 


»-f[ JEt+ r= « 

Since a + /3 and a/3 are both positive, a and £ must be both real and 
positive, or else conjugate imaginaries with positive real parts. In either 
case y tends to KE as a limit, as we should expect. 

W^e notice incidentally that if the circuit contained no capacity or 
self-induction the differential equation would be simply 


Ey = E. (5) 

Hence if the current has been found for simple resistances, self-induction 

and capacity can be allowed for by writing Lp + R + ~ f or R. For this 

Kp 

reason this expression is sometimes called a ‘resistance operator,’ and 
the method generally the ‘method of resistance operators.’ 


3.2. The Wheatstone Bridge method of determining Self-Induction . 
In this method the unknown inductance is placed in the first arm of a 
Wheatstone bridge; the fourth arm is shunted, a known capacity boinu 
placed m the shunt*. J 

First consider the ordinary Wheatstone bridge, the resistances of the 
arms being j?,, R„, R», R t ; let x be the current in R lt y that in 74, q that 
through the galvanometer, and b the resistance of the battery and leads. 



* Bromwich, Phil. Mag., 87 (1919), 407-419. Further references are given. 
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Then RiX-Rzy + Gg-Q, (1) 

R$x — R 4 y — ( R$ + R 4 + Gf) g = 0, (2) 

b (x + y) + -Ka# + R* {y + ff) = E, (3) 

and on solving we find 


9 


® + y 


00 


jffio j® 8 — Mi jfi>4 ^ + .fig + jSs + jf? 4 ) + (M l + jBq) (jS 8 + M^) 

If b is large compared with i? 2 and R 4 , we have nearly 

x+y^E/b, (5) 

/ 3 

and 


6? (jRl + J?2 + jBs + A> 4 ) + (JSl + M 2 ) (*Bs + J?0 


( 6 ) 


According to the result of 3.1, we can allow for the self-induction in 
the first arm by replacing R x by Lp+ R x . Let the arrangement in the 



Fig. 3. 


fourth arm be as shown (Pig. 3). The resistance of the main wire is R 4 , 
that of the shunted portion of it r. Suppose the shunt to have a 
resistance J3. Then the effective resistance of the whole arm would be 


R 4 -r + 


rS 
r + S 





( 7 ) 


If the shunt contains a capacity K ] we must replace S by 8 + 1/Kp. 
Hence in the formula for g we must replace R x by Lp + R x and R 4 by 

R zEe 

(r + ST) Kp + 1 * 

The result expresses the current through the galvanometer when the 
battery circuit is suddenly closed. 

It can be shown that in actual conditions g cannot vanish for all 
values of the time. A sufficient condition for this would be that the 
operator R 2 Rz-R 1 R 4 should be identically zero; then g would be 
identically zero whatever the remaining factor might represent. But 
with our modifications this factor becomes 

® 


3*2 
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Multiplying up and equating coefficients of powers of p to zero, wo 
find 

B i (r + S) = r‘, 

- LB, + (M-M) (r + 8) K+ R^K-O. (10) 

BaBs — B\B, — 0 (11) 

From the construction of the apparatus r ^ B t , r + 8. Thus (9) can 
hold only if r = B, and 8= 0. The shunt wire must be attached to the 
ends of B, and must have zero resistance. Substituting in (10) wo find 

L=B 1 B i X, (12) 

together with the usual condition (11) for permanent balance of the 
bridge. 

Actually these conditions for complete balance cannot be completely 
satisfied ; but for the determination of L it is not necessary that they 
should. Suppose the galvanometer is a ballistic one, and that it is so 
adjusted that there are no permanent current and no throw on closing 
the circuit. Thus 

Lim g = 0; f gdt= 0. , (13) 

t-*»oo J 0 


.(14) 


If g is expressed in the form f(j>)IF{p), we have 

F(S>) *F'(a) e * 

where all the a s, in the conditions of the problem, will have negative 
real parts. Then the condition that g shall tend to zero gives 

/(0) = 0 (ifi) 


' A - ls0 /„ 9 dt ~ ^ Jf' ( a ) ^ 

Equation (15) shows that the operational form of g contains p as a 
factor. Also we can write 


•g zoo 

o.F‘ (&)p — a ’ 


( 17 ) 


so that (16) is the limit of g/p when p tends to zero. The vanishing 

of Lim g and / gdt implies therefore that the operational form of <j 
J 0 

contains p 2 as a factor. Hence the modified form of - Il x R I(L must 
contain p* as a factor. Then (10) and (11) still hold; (9) no longer 


holds. We now have 

R 2 Rs-R x R 4 = 0, (18) 

LRt=R x **K, ( 19 ) 


which gives the required rule for finding L. 
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3.21. Discharge of a condenser in one of two mutually influencing 
circuits. Suppose that we have two similar circuits, each with self- 
induction L and capacity K, but negligible resistance, and that the 
condenser in one carries a charge x 0 . The coefficient of mutual induction 
is M. The first circuit is closed; find the currents. 

If y denote the charges on the condensers in the two circuits, we 
have the differential equations 

Lx -J- + La?cc = 0, I ^ 

Mx + Ly + LoPy-0, J ^ ^ 

where ILL a 2 = 1 (2) 

Initially x = x 0i y=Q, x~y = 0 (3) 

The subsidiary equations are 

L(p* + **) x + Mp 2 y = Lp*x 0 , l m 

Mp\x + L(p* + a:*)y~Mp*x<>.] w 

Then solving by algebra 

JJ(p* + a»)~My LMaV () 

x ~jfip*+a i y-My p J ~ wy+ay-My y ; 
Put M = L/3. Then 

v y + a?y - my = d (y + a 3 - py) y + a-+py) ; ...(6) 

p 2 ( 1 — fif) + a 3 , 

^ ~ TCI - p) f + a 2 } {(1 + P)f + yP 
r a-^)p a , 1L tM.lv 

~ L(1 -P)ir + o: i {\ + p) 1 f + a 2 J “ 2 

(,) 

/3a 2 , 

y ~ {a - p)f + «*} {(i + p) y + « 3 } p Xo 

1 r i-/3 i + p -] . 

“ " 2 L(l-/S)j0 2 + a 2 (1 + P)p* + aU P x ° 

--*[*“ (8) 

W8 writ© 7 ^s“T + 8; (« 

the solutions take the forms 

x = x 0 cosytcos St 9 | 

y — Xq sin yt sin 5#. J ^ ' 
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Thus when /8/a is small the disturbance consists of a rapid oscillation 
with beats, the oscillation being transferred from one circuit to the other 
in each quarter-period of the beats. 


3.3. The Seismograph. In principle most seismographs are Euler 
pendulums — pendulums with supports rigidly attached to the earth, so 
that when the earth’s surface moves it displaces the point of support 
horizontally and disturbs the pendulum. The seismograph differs from 
the Euler pendulum as considered in text-books of dynamics in two 
ways : instead of being free to vibrate in a vertical plane, it is con- 
strained to swing, like a gate, about an axis nearly, but not quite, 
vertical, so that the period is much lengthened; and fluid viscosity or 
electromagnetic damping is introduced to give a frictional term pro- 
portional to the velocity. The displacement of the mass with regard to 
the earth then satisfies an equation of the form 

x + 2kx + n 2 x = X£, (1) 

where £ is the displacement of the ground, and k, n } and X are constants 
of the instrument*. Suppose first that the ground suddenly acquires 
a finite velocity, say unity. Then £ jumps from 0 to 1, and therefore cu 
from 0 to X. The initial conditions are therefore 


x = 0; £ = X, (2) 

and our subsidiary equation is 



( p 2 + 2kp + n 2 )x = \pH(t) 

( 3 ) 

Put 

p 2 + 2k p 4* n 2 = (p + a) (p + /?) 

( 4 ) 

Then 

-(I 

( 5 ) 


= 0 when t <0 

( 6 ) 

and 

= (jr^ t —er at ‘) when t>0 

( 7 ) 


The recorded displacement x therefore begins by increasing at a finite 

i 


rate X, reaches a maximum X 
tends asymptotically to zero. 



after a time 


log and then 


* Some instruments, such as that of Wiechert, are not on the prinoiple of the 
Euler pendulum, but nevertheless give an equation of this form. 
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If a and fi are real, and /*? less than a, we see that the behaviour after 
a long time depends mainly on e~ pt ; now as the experimental ideal is 
to confine the effects of a disturbance to as short an interval afterwards 
as possible, we see that we should make j3 as large as possible. But 


and for a given n, /3 is greatest when K-n. This is the condition for 
what is called aperiodicity ; the roots of the period equation are equal, 
real, and negative. Many seismographs are arranged so as to satisfy this 


condition. The solution is then 

® 

= 0 when t <0... (10) 

= Xte~ nt when £>0 (11) 


The maximum displacement is now at time 1/n after the start, and is 
equal to Xjen. 


If k < n, we can put 


(12) 

*?o 

1 

** 

II 

So 

Then the solution is 

%=^e~ Kt miyt y 

(13) 



Fig. 4. Recovery of seismographs with K~n and K—njJ 2 after same impulse. 
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and the motion does not die down so rapidly as for K = n. The aperiodic 
state therefore gives the least motion after a long interval for a given 
value of n . 

In practice, however, k is usually made rather less than n; in the 
Milne-Shaw machine, for instance, k is about 0’7 n. The motion is then 
oscillatory, but the ratio of the first swing to the second is e K7r/v , or about 
20. But x vanishes after an interval irjy from the start, or about 4/w, and 
ever afterwards is a small fraction of its first maximum. The reduced 
damping effect after a long time is considered less important than the 
quick and complete recovery after the first maximum. The time of the 
first maximum is 11 jn from the start, as against Ijn for the aperiodic 
instrument and 1‘57/w for the undamped one. 

3.31, The Galitzin seismograph is similarly arranged, but the motion 
of the pendulum generates by electromagnetic induction a current, 
which passes through a galvanometer. If x is the displacement of the 
pendulum, and y that of the galvanometer, the differential equations 
are 

X + 1k x X + UiX = ( 1 ) 

y + 2k 2 ?/ + n?y = fix (2) 


Supposing the ground again to start suddenly with unit velocity, we 
have 



Fig. 5. Recovery of Galitzin seismograph after impulse. 
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As a rule the two interacting systems are so arranged that k and n are 
the same for both, and both are aperiodic, so that * = n. Then 




(4) 

=v K!i e ~ n< ) when#>o 

(5) 



(6) 


The indicator therefore begins to move with a finite acceleration, instead 
of with a finite velocity as for the pendulum. The maximum displace- 
ment follows after time (3 - */3)/rc = l’27/rc, the mirror passes through 
the equilibrium position after time S/n, and there is a maximum dis- 
placement in the opposite direction after time 4‘73/rc. The mirror then 
returns asymptotically to the position of equilibrium. The ratio of the 

two maxima is e /( 2 + = In comparison with an instrument 

such as the Milne-Shaw, recording the displacement of the pendulum 
directly, the Galitzin machine gives the first maximum a little later, 
the first zero a little earlier, and the next maximum displacement is 
larger in comparison with the first maximum. 

In an actual earthquake the velocity of the ground is annulled by 
other waves arriving later ; the complete motion of the seismograph is 
a combination of those given by the separate displacements of the ground. 

3.4. Resonance. A simple pendulum, originally hanging in equilibrium, 
is disturbed for a finite time by a force varying harmonically in a period 
equal to the free period of the pendulum. Find the motion after the 
force is removed. 

The differential equation is 


x + n 2 w “/sin nt 

(1) 

=f n p ~ 

J p i + ?i a 

(2) 

The solution is then 

(3) 


nothing having to be added on account of the initial conditions. 
To evaluate this, we notice that 
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Differentiating with regard to n, we hare 

- ■ / 2\ 2 = - * COS nt - sin nt, (5) 

(jgt? 2 4- tzj n n? ’ w 

and as = (sin - nt cos raf) (6) 

Suppose the disturbance acts for time nr/n, where r is an integer. At 
the end of this time 

®=o (7) 

The subsequent motion is therefore given by 

X= ~ 2tf - eoa(«rt-rr) 

= -gfcOSM* (8) 


3.5. Three particles of masses §^?ra, and ?ra, in order, are attached 
to a light stretched string, the ends of the string being fixed. One of the 
particles of mass m is struck by a transverse impulse L Find the sub- 
sequent motion of the middle particle. (Intercollegiate Examination, 1923.) 

If x u x 2 , x s are the displacements of the three particles, P the tension, 
and l the distance between consecutive particles, we find in the usual 
way the three equations of motion 


#i = - A. (2#! - # a ), 

(1) 

fir &2 ^ — * A ( — ffii + 2x 2 #?a)j 

(2) 

#a = -M“#2 + 2tf 8 ), 

(3) 

p 

where A = — 

mi 

(4) 

Initially all the displacements are zero, 0, x x ~ I/m. Hence the 

subsidiary equations are 


(p 3 + 2A) a>! - Aa? a = p I/m, ... 

(5) 

— X#! + (|rj- p 2 + 2 A.) os 2 — Afl? s = 0, 

(6) 

- Aa? a + (jp 3 + 2A) co t = 0 

(7) 

Hence *, ^ +2i *.;®.-^, + - sii e i + y ,P 2X ^ 

(8) 

and (|£.p*+2A-- = Kp 1 . 

\20^ p 3 + 2\J j> a + 2A m 

(9) 

Multiplying by p 3 + 2A we find 


(7p* + 4A) (3p s + 10A) * a =20A^ I/m, . . 

(10) 
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20A.JO I 

X * ~ (7p* + 4X) (3jo a + 10X) to 

_i£_l 

7jo 2 + 4X ‘6p l> + 10XJ 


= 10 / 

29 m 

where a 2 = -£A; /? a = Xp-X (12) 

We notice that the mode of speed V(2A) is excited, but does not need 
to be evaluated because it does not affect the middle particle. 


f sm at _ sin p t' 

la 8 




.(ii) 



3.6. Radioactive Disintegration of Uranium . The uranium family 
of elements are such that an atom of any one of them, except the last, 
is capable of breaking up into an atom of the next and an atom of 
helium*. The helium atom undergoes no further change. The number 
of atoms of any element in a given specimen that break up in a short 
interval of time is proportional to the time interval and to the number 
of atoms of that element present. If then u, # 2 , ... x n are the numbers 
of atoms of the various elements present at time t , they will satisfy the 
differential equations 

du 1 

dt = ~ KU ' 


dx x 

^.-KU-K^u 

dx a } 

— K \^i — K a#aj 


(i) 


dXm 

dt ** Kn ~ l ■ 

j 

Suppose that initially only uranium is present. Thus when t = 0, w = w 0 , 
and all the other dependent variables are zero. Then the subsidiary 
equations are 

(p + k) u — pu a , ) 


(p + k,) x x = m, 

(p + K„) 0 a = K X X U 


00 


(.P + K n-l) ®n-i — K »-2*®»-a> 
px n ~ Xn—i • , 


We neglect /3-ray products, for reasons that will appear later. 
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The operational solutions are 

pih K pu o / 

U ~ p + K * ^ (jp + *) (#> + *i) * 

- - . 
a_ G> + «)Cp + 0(p+K»)’ 


(») 


71 ^ O + «) (p+*i) — (jp+ **-0' J 

These are directly adapted for interpretation by the partial-fraction 
rule. In fact 

u = u^er* 1 ; - e~ K \ l ) ; 

Kj — K 

$2 — KK^V/q ip T“7 ” 6 Kt + 5 -V"V v & 

l( K l “ K ) \ k 2 ~~ K ) ( K -« l ) (**-*0 /JX 

i 'i y •••(4) 




K i ■ 






Of all the decay constants k is much the smallest. If the time elapsed is 
long enough for all the exponential functions except e~ Kt to have become 
insignificant, these results reduce approximately to 


u=u 0 e Kt ; = -- u^er Kt \ x 2 = — u 0 e~ Ki ; ... a n ~u Q (1 -e~ Kt ). ...(5) 

K x K 2 

With the exception of the last, the quantities of the various elements 
decrease, retaining constant ratios to one another. 

On the other hand, if the time elapsed is so short that unity is still 
a first approximation to all the exponential functions, we can proceed 
by expanding the operators in descending powers of p and interpreting 
term by term. Hence we see that at first cc x will increase in proportion 
to t, #a to t\ and x n to 

In experimental work an intermediate condition often occurs. Some 
of the exponentials may become insignificant in the time occupied by 
an experiment, while others are still nearly unity. We have 

«>r = Kj ~T “ Kr ~ l [P~ l Mr-i - K r p~~X r _ x + . . .] ( 6 ) 

and if K r t is small we can neglect the second and later terms in comparison 
with the first. Thus in this case 

X r - K r _i p-'Xr-x (7) 



I 



i 


> 
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If a r -i is of the form t 8 , we can put 

_ _ _ K r— l £ j _ K r-lSl * * /n\ 

^ (i? + «r)i? 8 i? 8 ' 1 ’ 1 jP + *r ^ 8+1 ^ ' 

If K r t is small, we can replace the exponential by unity and confirm (7). 
If it is great, 

'p~^ e r K r ^ = J* e~ K r f dt = ™ + 0(^“V), 


and on continuing the integrations 

0~~K r t = /p~ 8 A — 

K r 


1 if 
K r S \ “ 


Thus = — # r -x (9) 

K r K r 

Classifying elements into long-lived and short-lived according as K r t is 
small or large for them, we find that the quantity of the first long-lived 
element after uranium is proportional to t, the next to t 2 , and so on. All 
/3-ray products are short-lived when t has ordinary values. 

Radium is the third degeneration product of uranium. In rock 
specimens the time elapsed since formation is usually such that the 
relations (5) have become established. As a matter of observation the 
numbers of atoms of radium and uranium are found to be in the constant 
ratio 3*58 x 10" 7 . This determines k/k q . Also the rate of break-up of 
radium is known directly : in fact 

1/k 3 = 2280 years. 

Hence 1/k = 6*37 x 10° years. 

This gives the rate of disintegration of uranium itself. 

A number of specimens of uranium compounds were carefully freed 
from radium by Soddy, and then kept for ten years. It was found that 
new radium was formed; the amount found varied as the square of the 
time. This would suggest that of the two elements between uranium 
and radium in the series one was long-lived and the other short-lived. 
Actually, however, it is known independently that both are long-lived. 
The first, however, is chemically inseparable from ordinary uranium, 
and therefore was present in the original specimens ; initially, instead 
of %i - 0, we have 

K 

K 1 

For the next element, ionium, we have 

— KU 0 t, 
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the variation of x x being inappreciable in the time involved. Also 

X S = Kg p" 1 #* = %KK 2 U 0 t*. 

Soddy* found that 3 kilograms of uranium in 10*15 years gave 
202 x 10“ 12 gm. of radium. Hence, allowing for the difference of atomic 
weights, 

a ? 3 /« 0 = 7 *1 * 10 ’ u 

and /f a = 8*64 x 10” fl /y ear 3 l/*«a= 1*16 x 10 B years. 

This gives the rate of degeneration of ionium. Soddy gets a slightly 
lower value for l/x a from more numerous data. 


3.7. Small Oscillations in Dynamics . Suppose that we have a 
dynamical system in equilibrium or steady motion, with a Routh’s 
modified Lagrangian function of the form 


2jK — a^sX r CPs -I" Gra&'r&'st (1) 

where x T are the departures of the non-ignorable coordinates from their 
values in the steady motion. The friction, if any, is represented by a 
dissipation function F given by 

2F=f ra x r x s (2) 

By symmetry there is no loss of generality in taking 

Mrs = O'er j &ra = Car j fr% —far (3) 

The equations of motion for small disturbances are 



dt \dx r ) dx r dx r 5 

( 4 ) 

that is 


( 5 ) 

where 

““ Qsr Qra + fra 

(6) 


If an external force S r is applied to the coordinate x r > the equations 
have the form 1.72(1). If . the system starts from rest and only one of 
the S r differs from zero, the subsidiary equations are 

e ma x a = 0 = (m = r) ( 7 ) 

Writing A for the determinant formed by the c’s, and E rs for the minor 
of e r8 in this determinant, we have the operational solution 

«.-!=«■ (8) 


If all the b n = 0, the determinant A is symmetrical, so that 

J&.-J&T 


* 


& 


* Phil. Mag. (6) 38 (1919), 483-488. 


( 9 ) 
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Then a given force S r applied to the coordinate x r will produce precisely 
the same variation in x 8 as the same force would produce in x r if it was 
applied to x 8 . Thus we have a reciprocity theorem applicable to all 
non-gyroscopic and frictionless systems. Friction does not affect the 
result if it is expressible by a dissipation function F, for f rs =/«* and 
the relation (9) still holds. 

If the forces reduce to an impulse at t = 0, so that S r can be replaced 
by pJ ri the solution becomes 

H , (io) 

We can evaluate the initial velocities by expanding in descending 
powers of p. The first term is 

= ( 11 ) 

where A is the determinant formed by the ct’s, and A rs the minor of 
a n in it. Hence the initial velocities, found by operating on this with 
p , are 



Now the constants a r8 are merely twice the coefficients in the kinetic 
energy, which is a quadratic form. Hence the determinant A is sym- 
metrical whether the system is gyroscopic or not, and the reciprocity 
theorem for impulses and the velocities produced by them is proved. 


3.71. The subsequent motion can be investigated by interpreting 
according to the partial-fraction rule. But let us consider the simple 
case where the system is non-gyroscopic and frictionless. Then 

A = J.n (^ 2 + a a ), (1) 

where the as are the speeds of the normal modes. Then 


x$ = X 


.. E»( r O 

a(dkfdp 2 ) P 'i=~ a ? 


J r sin at, 


,( 2 ) 


where E rs (- a?) denotes the result of putting -a 2 for p 2 in E r8 . The 
contribution of the a mode to the initial rate of change of cc 8 is therefore 


v a * 


■%.(-<*) r 

(dkjdp^p r 


,( 3 ) 


An immediate consequence of the presence of the factor E r8 (- a 3 ) in 
the numerator is that if an impulse is applied at a node of any normal 
mode, that mode will be absent from the motion generated. 
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Another illustration is provided by Lamb's discussion* of the waves 
generated in a semi-infinite homogeneous elastic solid by an internal 
disturbance. The normal modes of such a system include a type of 
waves known as Rayleigh waves. These may be of any length, and 
involve both compressional and distortional movement; if the depth is 

the amplitude of the compressional movement in a given wave is 
proportional to e - as , and that of the distortional movement to e~P* 9 
where a and ft depend only on the wave-length. Lamb found that if 
the original disturbance was an expansive one at a depth/, the ampli- 
tude of the motion at the surface contained a factor but if the 
original disturbance was purely distortional, the corresponding ampli- 
tude contained a factor These factors are the same as would occur 
in the compression and distortion respectively at depth/ in a Rayleigh 
wave with given amplitude at the surface. 

3.72. In the discussion of the oscillations of dynamical systems about 
equilibrium the ordinary method of seeking for solutions of the form 
® a = K0 yi meets with a difficulty when the determinantal equation for 
the periods has equal roots. This does not arise in the present method. 
If the system is not dissipative and the roots are unequal, we know that 
the zeros of the minor of any element in the leading diagonal separate 
those of the original determinant; hence in the limit, if the determinant 
has a factor (p 2 + a 2 ) r , the minor of any element in the leading diagonal 
has a factor (p 2 + a 2 /- 1 , and therefore every first minor contains this 
factor. Hence when we evaluate the operational solution for given 
initial conditions as in 1.72(7), in the absence of external forces, we 
shall have 

jgj 

and the factor (p ! + a 3 ) r ~ 1 cancels from the numerator and denominator 
of E r J A. Hence only a single factor (p* + a 3 ) remains in the denomi- 
nator, and the corresponding partial fraction contributes only trigono- 
metric terms cos at and sin at to the solution. 

In dissipative or gyroscopic systems, b„ is not zero, and the root- 
separation theorem may not hold. Then equal roots do not necessarily 
imply a corresponding factor in E rm , and the operational solution may ' 
still have a repeated factor in the denominator. Then terms like te~ ai , 
t cos at, t sin at will occur in the solution. A simple instance of this has 
already appeared in the aperiodic seismograph discussed in 3.3. 


* Phil. Trans., A, 208 (1904), 1-42. 
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CHAPTER IV 

WAVE MOTION IN ONE DIMENSION 

4 . 1 . In a large class of physical problems we meet with the differen- 
tial equation 

•••« 

•where t is the time, x the distance from a fixed point or a fixed plane, 
y the independent variable, and c a known velocity. Let us consider the 
solution of this equation first with regard to the transverse vibrations 
of a stretched string. In this case we know that 

c 2 = jP/«&, (2) 

where P is the tension and m the mass per unit length. Write p for 
d/dty and D for d/dx. Suppose that at time zero 

y=f {«=)■, |=^(4 (3) 

where /and F are known functions of x; that is, we are given the 
initial displacement and velocity of the string at all points of its length. 
Then we are led by our previous rules to consider the subsidiary 
equation 

( p* — &D' 1 ) y -£>*/(^) +P F(x) (4) 

or y ~ W + jjrz jy* ^ 0 27 ) (®) 

But P. + _JL\ 

%r - c l D u z \p-cD p + cDJ 

= + ( 6 ) 

P _ ( P . _ P \ 

p~ - rJP 2 cJJ \p - cD p + cl)/ 

= ( 7 ) 

In (7) the l/D has been put last in consequence ot our rule that opera- 
tions involving negative powers of D must be carried out before those 
involving positive powers. Hence 

m = i (« em + 

- J {/(® + Gi ) +f - C 0}: (3) 

J 4 
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1 p fX+Ct rX—Ot -I 

= ¥ c[L * v >**-1 

1 {X -f- Ct 

<•> 

and therefore 

1 rX+at 

2'=il/(« f + ^)+/(«-^)} + 9r/ F{x)dx. (10) 

*CJx-Ct 

This is D'Alembert’s well-known solution. 

This cannot, however, be the complete solution. Equation (1) has 
been assumed to hold at all points of the string; but if any external 
forces act these must be included on the right of the equation of motion. 
In such problems the ends of the string are usually fixed, and reactions 
at the ends are required to maintain this state; in the complete equa- 
tions of motion these reactions should appear. They are unknown 
functions of t, and therefore necessitate a change in the mode of solution. 
This corresponds to a real indefiniteness in the solution itself; for if the 
string stretches from #=0 to x +ct) is unspecified if x + ct>l, and 

f(x~ct) is unspecified if x~ct is negative. The solution is therefore 
incapable of application to a specific problem unless we can obtain 
a means of specifying /(#) and Fix) outside the original range. 

But D’Alembert's solution can be adapted so as to fulfil all the con- 
ditions. We notice that the solution consists of two waves travelling in 
opposite directions with velocity c. If the initial disturbance is confined 
to a region within the string separated by finite intervals from both 
ends, it will take a finite time before D’Alembert's solution gives a 
displacement at either end ; hence no force is required to maintain the 
boundary conditions, and the solution will hold accurately until one of 
the waves reaches one end. 

Again, the initial disturbance is specified only for points within the 
length of the string, that is, for values of x between 0 and say; and 
by the last paragraph D’Alembert's solution would be complete if the 
length were infinite. If then we consider an infinite string stretching 
from - oo to + co , and disturbed initially so that fix) and Fix) are both 
antisymmetrical about both x—0 and x~l,y will vanish ever after at 
these points, so that no force is required to keep them fixed, and 
D’Alembert’s solution for such an infinite string will therefore be 
equivalent to that for the actual string from x - 0 to x = l. Thus we have 
a formal rule for finding the position of the string at any subsequent 
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time : consider an infinite string with the actual values of the initial 
velocity and displacement between # = 0 and x = /, but with the dis- 
placement and velocity outside this stretch so specified as to be 
antisymmetrical with regard to both % = 0 and x~l\ then D’Alembert’s 
* solution for the infinite string will .be correct for the actual string for 
the same values of x. 


4.2. We may approach the problem in another way. Talcing the same 
subsidiary equation 

W 


except that the effects of the initial velocity will not be considered at 
present, let us solve with regard to x as if p were a constant. The 
justification of this procedure is that 9/9# is commutative with regard 
to both differentiation and integration with regard to t. The boundary 
conditions are that y = 0 when # = 0 and when x = l. Using the method 
of variation of parameters, we assume that the solution is 


y = A cosh^— + B sinh , (2) 

where A and B are functions of x subject to 


A ' cosh — + B f sinh^ 7 = 0 (3) 

go w 

Substituting in the equation (1) we find 

A ' sinh^’ + B' cosli^-- = -£/(*) (4) 

c 0 c v y x ' 

Hence 

A ' = ~f(x) sinh — ; B' = — cosh ~~ (5) 


Now y must vanish for all time when x = 0 ; hence A (0) =0. Thus 



A ~ P l f(0 sinh — dt 

Jo c c 

(6) 

Also y must vanish when x= L Hence 



B(l) sinh $ = — A ( l) cosh ~ 

(7) 

and 

B(l)=~ coth f -/ (£) sinh — d£, — 

G Jo G C 

(8) 

giving 

B(tv) = - coth — f -/(£) sinh^- d£ 

C Jo c c 



+ I f/(£)oo8b*fdf. 

Jb C o 

(9) 
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In all we find 


y = jT|/(0 sinh^(coshf - coth^ smh fj % 

+ I f/(€ sin]l ~~ (cosi-^ -coth-^smh-^W 

JxC C \ G C CJ 

L 


* ™ 

Before carrying out the integration with regard to £ we can interpret 
the integrand by the partial-fraction rule. Bach integral vanishes when 
p = 0, and therefore no constant term appears in the solution. But 
&mhpllc vanishes when 

p1 z 


= rur. 

c 

where r is any integer, positive or negative. Hence 

nx . . ru rf . , l -a 

* . i /7 \ / - 7 -smn— r sinhn7r — y- 

/ (/-aQ/c _ l l l_ 

c c smhpl/c r=—oo crt-n- , 


.(id 


fvrrotll 


-j- cosh (nir) - 

6 C 


® 2 . rirf . nrc# 

= 2 rsm -j- sm -y- cos — j- , 

r=l Lit L 


.( 12 ) 


By symmetry the corresponding factor in the second integral has the 
same value. Hence 


rivet 


2 2 [ l .... . nr£ . . ritx 

S y J /(f) sm -y if . sm y cos 

iution given by the method of normal c 
ain the Fourier sine series 

/(•) = 5 J jf ‘ /(|) sin sin ™ d£ (14) 


(13) 

This is the solution given by the method of normal coordinates. Putting 
t zero, we obtain the Fourier sine series 


The effect of an initial velocity can be obtained by a similar process. 

But in practice the solution by trigonometrical series is not often the 
most convenient form unless we are satisfied with a knowledge of the 
free periods. It usually converges slowly; but what is more serious 
is that its form suggests little about the nature of the actual motion 
beyond the fact that it is periodic in time 2 Ijc. To find the actual 
form of the string at any instant it is necessary to find some way of 
summing the series, which may be rather difficult. A more convenient 
method is often the following*. We have seen that the interpretation of 
* Heaviside, Electromagnetic Theory , 2, 108-114. 
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any operator valid for positive values of the argument is equivalent to 
an integral in the complex plane, the path of integration being on the 
positive side of the imaginary axis. Then a factor cosech pl/c in the 
operator can be written 

and when interpreted by Bromwich’s rule this will give rise to a series 
in powers of e~ 2xl!o in the integrand. But since the real part of A is 
everywhere positive on the path, the series is absolutely convergent, 
and integration term by term is justified. Hence it is legitimate to ex- 
pand the operator in this way and to interpret term by term. Then for 

amhp£/c8mhp(l-ai)/c = } „ _ ?)/0 . _ ^ 

amh. pl/c 1 K ‘ 

(1 — -“)/<=) (l+ e -W0 + e -4j>?/0+ ...) ; ...(16) 

and for i>£, 

sinhjP^/c ainhy (l- £)/c = , _ x)/o ( _ e _ 2px/c) 

sinh pl/c 2 v ' 

(l_g-2p(J-f)/0) (1 + g-ipl/d + e -wic + ...). ...(17) 

On multiplying out either expression we obtain a series of negative 
exponentials of the forms e~* (a ~ f)/<! , with a > ij, or (?-»)/«, with £> a. 

Then j*£f®e-*0-Wdi = jf(£)de~^-» (18) 

But according to our rule e-v ( a ~M c is zero when t~(a — £)/c is negative 
and unity when this quantity is positive. Hence the integral (18) is 
zero unless a — ct lies between £ x and £ 2 , and then it is equal to f(a-~ ct ). 

Similarly f^M) r**—* * — Jf(() de -*«-■* (19) 

and = 0 if $ > ct + a 

= 1 if £^ ct + a } 

so that (19) is zero unless ct + a lies between £ x and £ a , and then is equal 
to f(ct + a). 

It follows that at any instant the effect of the initial displacement 
/(£) at the point £ is zero except at a special set of points where one 
or other of the quantities a ± ct is equal to £. Since a does not involve 
£ } we see that this way of expressing the solution reduces it to a set of 
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waves moving in each direction with velocity c. The first three factors 
in (16) give 

| { e -p(p ~£)/e _ e -p(x+£)lc _ g-#{(aJ-£)+2 ( l-<e)}/c + g-#{(aH-£)+2 (X— 

( 20 ) 

for values of x greater than £ The first term gives J/(£) at time 
(®-£)/c 3 the second -A/(£) at time (#+£)/<?, the third -&/(£) at 
time {# - £ + 2 ( l-x)}/c and the fourth J/(£) at time \x +■ £ + 2 (l-x)}/c* 
The first term represents the direct wave from £ to x } the second the 
wave reflected at x = 0, the third that reflected at x = <?, and the fourtli 
one reflected first at x - 0 and then again at x = l. The term rW in 
the last factor of (16) will give four further and similar pulses, each 
later by 2 l/c than the corresponding pulse just found. These are pulses 
that have travelled twice more along the string, having been reflected 
once more at each end. Similarly we can proceed to the interpretation 
of later terms as pulses that have undergone still more reflexions. 

The part of the solution arising from values of £ greater than x may 
be treated similarly. The interpretation in terms of waves is exactly 
the same. 


4.3. As a particular case let us consider a string of length l , origin- 
ally drawn aside a distance y at the point x = b, so that initially it lies 
in two straight pieces, and then released. Then 

Jfo = yx/b O^x^b) 

yQ-y(l~x)l(l-~b)- 

and the subsidiary equation is 

&y p* p 2 

da? oo 

The solution that vanishes at the ends is 

y = y g + A sinh-2— sinh £(l- b) 


0 ^ x ^ b j 


+ j 4 sinh-^ sinh b^x^l j 


•GO 


'b 

l~X 

c c 

where the constants have been chosen so as to make y continuous at 
<*-b Also a discontinuity in dy/dx at this point would imply an infinite 
acceleration which cannot persist. Hence we add the condition that 
%/3# shall be continuous, which gives 

’ (j + ih) * *7 {“‘ h f ** f (*- *0 4 sinh^ oonli f (J-4)} „ o 

« 
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and on simplifying 

A = - - 7 ,/ 7) cosech 

p b (l— b) o 

Thus 

y = f ^ g _ £ a mh^/csinhp (l-b)jcx Q ^ x%b 

J b(l-b){ l p smh pl/c J 

^ b(l-b) \l K ' p smhjpZ/c J 


..(5) 

..( 6 ) 


Interpreting by the partial-fraction rule, we notice that the contribu- 
tion from p = 0 just cancels the term independent of p, while the rest 
gives* 




2 r)P ” 1 . 

~ “t —t 


0 ^ < £ 


( 7 ) 


Alternatively, we can interpret the solution for 0 < # < b in exponentials, 
for 

sinli pxjc sinh jo (l—b) /c = ^_ p (s _ 8)/c {1 _ 
sinhjpZ/c - v 

„ e -2p(l-1>)/c') (l + 0 -S^/c + 0 -4pZ/c + .(8) 


and ~=c£ £ > 0 (9) 

P 

Then 0-P(b-x)io^ ct -^-ai), (10) 

P 

when the right side is positive, and otherwise is zero. Thus y retains 
its initial value »?#// unaltered until time (!) -x)/c, when the wave from 
the part where x > b begins to arrive. After that we shall have 

» - WH) fr * - i c* - *♦ •>} ' STPT) ((I - 7) i < c ‘- «} • 

( 11 ) 

This solution remains cori’ect until ct = b + x, when the wave reflected 
at x — 0 begins to arrive. Thenceforward 


Hence 


(<r* >(*-»)/« - e -p { <’+»)/”) _ = [ct - Q) - oi)) ■ 

l — h 


■{ct-(b + ai)} = 2x. ...(12) 


y 


v j 

b(l-b) 


fl-h \ 

("7”®“®/ = - 


•tjX 

i-y 


.(13) 


The part of the string reached by the first reflected wave is therefore 
parallel to the original position of the part where b<x<l. 


Bayleigb, Theory of Sound , 1, 1894, 185. 
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When ct = (b - as) + 2 (l - b'j, a wave reflected at x -l arrives, and 

y = b (1- b ) {(2 ~j) a, + H ct + b )~ 1 } (14) 

This holds until ct = (b -f cc) + 2 (l — b) ; in the next phase 





Fig. 6. 


When ct - 21, the whole of the string is back in its original position; 
the term in rW then begins to affect the motion, and the whole 
process repeats itself. We see that at any instant the string is in three 
straight pieces. The two end pieces are parallel to the two portions of 
the string in its initial position, and are at rest. For the middle portion 
the gradient dy/dw is the mean of those for the end portions, and the 
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1 Yfv(j • mm m 

transverse velocity is ± - ^ jfzft j ■ The middle portion is always either 
extending or withdrawing at each end with velocity c . (Pig. 6.) 


4„4- A uniform heavy string of length 2 1 is fixed at the ends. A 
particle of mass m is attached to the middle of the string. Initially the 
string is straight. An impulse J is given to the particle. Find the sub- 
sequent motion of the particle. (Of. Rayleigh, Theory of Sound, 1, 1894, 
204.) 

The differential equation of the motion of the string is 


dp W 


- 0 . 


( 1 ) 


Take as zero at the middle of the string. By symmetry we need consider 
only the range of values 0 ^ as < l Suppose that the displacement of the 
particle is rj. When t = 0, y and dy/dt are zero except at as - 0. The 
subsidiary equation therefore needs no additional terms. The solution 
is therefore 


y = v 


sinh ff ( l—as)lc 
sinli pi jo 


( 2 ) 


If now P be the tension in the string, the equation of motion of the 
particle is 



If p be the mass of the string per unit length, 

P = (4) 


On account of the initial conditions the subsidiary equation for the 
particle is 



m * fri= ‘ 2P \£L 0 +Jp 



~ - S Prj £ COtll — + Jp 

c c 

(5) 

Hence 

Jo 

mpe + 2 P oo tli pl/c 

(6) 

Put 

2 pi ^ mass of the string 
m mass of the particle * 

(7) 

Then 

b3 

li 

(8) 
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To interpret this solution by the partial-fraction rule, we -recollect that 
the system is a stable one without dissipation, and therefore all the zeros 
of the denominator are purely imaginary. Putting 

pljc = i o) (10) 

we see that the zeros are the roots of 


<o = A COt (o (11) 

There is a root between every two consecutive multiples of t, positive 
or negative; and the roots occur in pairs, each pair being numerically 
equal but opposite in sign. Then 


U l 

me Lwc(l + k cosec 2 <o) (Ijc) 




.( 12 ) 


the summation extending over all positive and negative values of o>, 


9 77 I 

= — 2 — 7 r— T — n sin (<*)c£/0, (13) 

me o) (1 + k cosec a>) v J v ' 

the summation being now over positive values of <*>. 

If a root of (11) is Mr + X, where n is great and X <-*, we have 

(wtt + X) tanX = k (14) 

and \=k/mr (15) 


approximately. The series (18) converges like 2 - sin 2 a>, or 2 . It 

therefore converges fairly rapidly at the beginning, but more slowly 
later. Also the w’s are incommensurable, and therefore the motion does 
not repeat itself after a finite time ; thus the labour of computation 
would be great. 

To express the solution in terms of waves, it is convenient to change 
the unit of time to l/c , the time a wave takes to traverse half the length 
of the string. Also Jjm may be replaced by V. Then 

V V(l- e~ 2p ) 

^ p + k coth p (p + k) - {p - k) e~ 2p 


= V — 
P 

V 

p + k 




H 


2k 


2 k ( p - k) 
(j> + kf 1 


-*—•) 06) 


The first term is zero for negative values of the time; for positive values 
it is equal to 


V 

k 


(I-*-*). 


(17) 
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After time 2 the second term no longer vanishes. We have 

k _ 1 p p_ 

(p + k?~ k Jc(p + A) ( p + A) 2 

(18) 

and the second term is equivalent to 

- {1 - - A (# - 2) for *> 2 (19) 

The third term is zero for t < 4 ; for greater values it is easily found to 
be 

OF 

~ [1 - {1 + A (* - 4) + A 2 (t - 4) 2 } e-«W] (20) 

The process may be extended to determine the motion up to any instant 
desired. The entry of a new term into the solution corresponds to the 
arrival of a new wave reflected at the ends. 


4.5. A uniform heavy bar is hanging vertically from one end, and a 
mass m is suddenly attached at the lower end. Find how the tension at 
the upper end varies with the time. (Love, Elasticity, § 283.) 

If x be the distance from the upper end, and y the longitudinal dis- 
placement, y satisfies the differential equation 




a) 


where p is the density, E Young's modulus, and F the external force 
per unit volume, in this case pg . Put E/p = c 2 , and let the displacement 
of a particle before the weight is attached be y 0 - Then 



00 


When x = 0, y Q = 0 ; and if the bar be of length l , dy Q /dx = 0 when x-L 

jl® i 


Hence 


?/o 


_glx( , x\ 


.(3) 


After the weight is attached, we still have 

W dw‘ 9 to 3 ’ W 

while when £ = 0, y = y 0 and dy/dt^O. Hence the subsidiary equation is 


w 
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and the solution that vanishes with os is given by 

y-y§ — A sinh pas) c, (6) 

where A is independent of os. 

If w be the cross section of the bar, the equation of motion of the 
mass m is 


m w = mg - Em b 


CO 


the derivatives being evaluated at x = /. The subsidiary equation is 


mp*y = mg + mp*y Q -EmA } ( 8 ) 

which gives on substituting for y from (6) 

( p 2 sinh — + cosh-^) A = g (9) 

V c me c) * K J 

The tension at the upper end is 

gpl + EpAjc gpl + ^ qq^ (pl/c) + mep sinh (pl/c ) * 

If k be the ratio of the mass of the weight to that of the bar, 


1c-— — -id 

purl’ B*~* e* 

and the tension is 


( 11 ) 



We see that gmj t mk is the tension due to the weight of the bar alone, 
and gm/vr is the statical tension due to the added load alone. To evaluate 
the actual tension, we expand the operator in (12) in powers of e -$ l l c % 
Taking l/c for the new unit of time, we have 


2e~ p 


kp sinh p + cosh p ( kp + 1) - {kp - 1) e^ } 

kp + lL kp + 1 \kp + l) 


( 13 ) 

The first term vanishes up to time unity, and afterwards is equal to 

2 ( 1 -«-(*-«/*) ( 14 ) 

This increases steadily up to time 3, when the next term enters. Again, 


kp-1 __ . kp 2 kp 

(kp + iy l+ kp + i (kp + iy 

= - 1 + erW* + 2 (tjk) e-W 


( 15 ) 
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and the first two terms, when t> 3, are equivalent to 

2 (1 - <r «-«/*) - 2 [1 - -2{(t- 3)/fy *-(*- »)/*] 

= 2«-(«-s)/* [1 + ( 2/k ) (# - 3) - e-^] (16) 

This has a maximum when 

1 + 0-w = H (# - 3) (17) 

Equation (17) has a root less than 5 if 

4/*>l + 0-**, (18) 


which is an equality if h = 2*7. Thus for A = 1 or 2 the maximum tension 
will occur before t = 5. If A = 1 the maximum tension is when t = 3*568, 
and is equal to 3*266 gm/w, that is, 1*633 times the statical tension. 
If h = 2 the corresponding results are t- 4*368, 2‘520 gm/m, and 1*680 
times the statical tension. 

The third term enters at time 5, and afterwards is equal to 
2 [1 - fl-<* - 6 >/ fe - 2 {(* - 5 )/£} 2 0 -<*- 6 >/*]. 

If h- 4, the maximum stress is when t = 6*183, and is equal to 2*29 gm/m. 
The statical tension is 1*25 gm/nr, so that the ratio is 1*83. 

Love proceeds by a method of continuation, but the present method 
is much more direct, and probably less troublesome in application. 

4.6. To prove that the results given by the operational method, when 
applied to the vibrations of continuous systems, are actually correct, it 
is necessary to show that the solution actually satisfies the differential 
equation and that it gives the correct initial values of the displacement 
and the velocity at all points. The proof that it satisfies the initial 
conditions is difficult in the most general case, but Bromwich has gone 
far towards one. We have seen that in one of the simplest problems, 
that of the uniform string with the ends fixed, the verification that 
the solution is valid for the most general initial displacement is 
equivalent to Fourier's theorem. But for more specific problems the 
operational solution is equivalent to a single integral, and the direct 
verification that it satisfies the initial conditions is usually fairly easy. 
To show that it satisfies the differential equation, it would be natural 
to differentiate under the integral sign and substitute in the equation. 
But in practice it is usually found that the integrand, near the ends 
of the imaginary axis, is only small like some low power of 1/A. (the 
second in the problems of 4.3, 4.4, and 4.5), and consequently the 
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integrals found by differentiating twice under the integral sign do not 
converge. But we can proceed as follows. If a part of our solution is 

« 

and this integral is intelligible for all values of x and t within, a certain 
range, we have 


- Lim ^ t+h ~ ^ + 

df " iSf A a 


=1 £stJ/w ex p x (* + 9 


^ h -2 + e~ Kh dK 


Also 


c 2 ^ = c 2 Lim ~ ~ l ~ 

fo 2 a-o &h‘ 




a?\e w -2 +<r A ' l <iX 


/r 


X ‘ 


.( 2 ) 


.(3) 


The two integrals are identical before we proceed to the limit, and 
therefore their limiting values are the same. Hence 

0£ 2 0 da?’ 60 

the /“ 6<1Uati0n is satisfied - This ar gument does not 

tX“S; v “ * " “ presdUe » 

The argument breaks down at points where the second derivatives do 
not exist; as for instance in 4.3 at the points where the slope of the 
curve formed by the string suddenly changes. At these points there is 
• , C ? ntimuty m the ^ansverse component of the tension, so that the 
pomt has momentarily an infinite acceleration. This is why the velocity 

dven sLtrTr°T y W f 6n a WaV6 arrivea Tlle momentum of l 
J T. f n\ f J he atnng ’ ll0wever > varies continuously; the difficulty 

™« rf ^ repr “““ 0n b? * ***** -I"*** »ot of J 
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CHAPTER Y 

CONDUCTION OP HEAT IN ONE DIMENSION 

5.1. The rate of transmission of heat across a surface by conduction 
is equal to — kd Vfin per unit area, where V is the temperature, h a con- 
stant of the material called the thermal conductivity, and dn an element 
of the normal to the surface. Hence we can show easily that in a uniform 
solid the rate of flow of heat into an element of volume dccdydz is 

Jc V s V . dxdydz . 

But the quantity of heat required to produce a rise of temperature dV 
in unit mass is cdV ’ where c is the specific heat, and therefore that 
required to produce a rise dV in unit volume is pcdV s where p is the 
density. Hence the equation of heat conduction is 


^{pcV^-kV-V ( 1 ) 

If we put 

i/pc = A 2 , (2) 

W is called the thermometric conductivity, and the equation becomes 

(3) 


In addition, there may be some other source of heat. If this would 
suffice to raise the temperature by P degrees per second if it stayed 
where it was generated, a term P must be added to the right of (3). 

It is usually convenient to write p for 3/3 1, and 

p = h\f (4) 

The operational solutions are then functions of q ; but q must be expressed 
again in terms of p before interpreting. 

5.2. Consider first a uniform rod, with its sides thermally insulated, 
and initially at temperature 8. At time zero the end x = 0 is cooled to 
temperature zero, and afterwards maintained at that temperature. The 
end x = Z is kept at temperature 8. Find the variation of temperature at 
other points of the rod. 

The problem being one-dimensional, the equation of heat conduction is 

dv , 2 3 a r 
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while at time 0 , V is equal to 8. Hence the subsidiary equation is 


d 2 V 

pr~^=i>s, 

(2) 

» 

(3) 

The end conditions are that 

V = 0 when x = 0, | 

V -8 when x = L J 

(4) 

E« 

\ smn ql J 

(5) 

The integrand is an even function of q } and therefore 
function of p. It has poles where 

a single-valued 

ql-± t-rnr ; that is, p-- k 2 ?i 2 7i 2 /l 2 J . . 

(6) 


where n is any integer. But the negative values of q give the same 
values of p as the positive values, and therefore when we apply the 
partial-fraction rule we need consider only the positive values. The part 


arising from p = 0 is 

to 

The general term is 

_ ft sin hnwg-g)/; 

(- cosh (un) ( P/2#m ) 

= /S' — sin ^ e -wmm ( 8 ) 

iT # Mr l w 

and the complete solution is 

F= # ff + S — sin ^ 1 (9) 

L l n^lllTT l J ^ J 


If Trht^/1 is moderate, this series converges rapidly, and no more con- 
venient solution could be desired. It evidently tends in the limit to the 
steady value Sx/l. 

But if Trht^/l is small the convergence will be slow. In this case we 
may adopt a form of the expansion method applied to waves*. We can 
write (5) in the form 

F= 8 [1 - e-a® (1 - &-W-*)) (1 + + *r 4 « z +...)] (10) 

For if we interpret this as an integral along the path L, the argument 
of q is between ± \tt at all points of the path, and the series converges 

* Heaviside, Electromagnetic Theory , 2, 69-79, 287-8. 
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uniformly. Integration term by term is therefore justifiable, and we may 
interpret term by term. Now 

qx = xjp^jh 3 (11) 


and by 2.4 (16) 
Hence 


e -w = l — erf 


x 


2k$' 


.( 12 ) 


F=/Sf"erf (l - erf — ^ -(l - erf — ^ )~...1....(13) 

L 2 kfi \ 2 h£' \ 2hi J 


When w is great, 1 - erf w is small compared with e~ w \ If then x/2 

is moderate, but l/2h$ large, this series is rapidly convergent, and can 
in most cases be reduced to its first term. This solution is therefore 
convenient in those cases where (9) is not. 


5.3. One-dimensional flow of heat in a region infinite in both direc- 
tions . First suppose that at time 0 the distribution of temperature is 
given by 

( 1 ) 

We have just seen that the function expressed in operational form by 


e~ qx = 1 - erf 
satisfies the differential equation 

2 = Oj 


X 


dy 

dt 


dx 2 


( 2 ) 

(3) 


for positive values of x; and by symmetry it will also satisfy it for 
negative values of x , since the function and its derivatives with regard 
to x are continuous when # = 0. Also when t tends to zero this function 
tends to zero for all positive values of x, and to 2 for all negative values. 
It follows from these facts that the function 


K2-0=|[l+erf 2 ^] (4) 

satisfies the differential equation for all values of x and all positive 
values of t; and when t tends to zero the function tends to zero for 
negative values of x } and to unity for positive values; that is, to If(x). 
Hence this function is the solution required. 

Suppose now that the initial distribution of temperature is 


V=f(fib ( 5 ) . 

J 5 
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where f(x) is any known function. Then this is equivalent to 

F=f >)} (6) 

Ji- -co 

But the solution for positive values of the time that reduces to H (£-#) 
when ^ = 0 is 



Hence, by the principle of the superposibility of solutions, the solution 


of the more general problem is 

® 

® 

Put now £ = # + (9) 

Then V^f_J{jc + 2ht^k)e-^dk, (10) 

This is the general solution obtained by Fourier. 


5„4. If the temperature is kept constant at #=0, but the initial 
temperature is f(x) for positive values of x , we may proceed as follows. 
If we consider instead a system infinite in both directions, but with the 
initial temperature specified for negative values of x so that 

/(-*)«-/(*), (1) 

we see that the temperature at x = 0 will be zero at all later instants and 
the solution of this problem will fit the actual one. Hence 


® 

= i f ,/(» + 2fe*A) <r^A - 4= r /(- a: + 2/rf*A) e^’rfA 

V 7T ' - aj/2/it* v 7T «'a;/27iis 

( 3 ) 

If in particular the initial temperature is everywhere S, f(£) = S, and. 


i .x/mi 

V=S. -i f e- 


e-^dk 


- 8 erf -® . 

2M 


( 4 ) 
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Thus the solution of 5.2 is regenerated. In Kelvin's solution of the 
problem of the cooling of the earth, 5.3 (10) was adapted to a seini- 
infinite region in this way. 

The value of dVfix at the end # = 0 is obtainable by differentiating 
the solution valid for a semi-infinite region. In Kelvin's problem 

F= $ (1 - <r a ®), 



The same result is found by differentiating (4). 

We notice the curious fact that although the original exact solution 
for a finite region in 5.2 (5) is a single-valued function of jt ?, a square root 
of t appears in the approximate solution for a greatly extended region. 
The reason can be seen by differentiating the exact solution. It gives 

GfL-* 00 ** (6) 

which is a single-valued function of p . But when we use Bromwich's 
interpretation we find that 



which is again single valued; but if A^ is specified to be real and positive 
when A. is real and positive, it has a positive real part at all points on 
L , and therefore if l is great coth A %l/h is practically unity. The integral 
is therefore equivalent to 



which is our interpretation of p^S/h, and is equal to SlhJ(Trt). We could 
have started by specifying the sign of A^ to be negative when A is positive, 
but then A^ and coth A Hfh would both have been simply reversed in sign, 
and the same answer would have been obtained. 

5.5. Imperfect cooling at the free end of a one-dimensional region. 
With the initial conditions of 5.2, let us suppose that the end x =* Hs 
maintained at temperature 8 as before, but that the end x = 0 is not 
effectively cooled to temperature zero. Instead we suppose that it radiates 
away heat at a rate proportional to its temperature. At the same time 
heat is conducted to the end at a rate kd V/dx per unit area. These 
effects must balance if the temperature at the surface is to vary 

5 -* 
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continuously, so that instead of having F= 0 at the end as before we 
shall have a relation of the form 


dV 

fa? 


aV = 0 at x- 0. 


a) 


The equation 5.2 (3) is unaltered, and the operational solution is 


F=tf{l-^sinhg (*-<*)} (2) 

where A is to be determined to satisfy (1). Hence 

qA cosh ql - a (1 - A sinh qt ) = 0 (3) 

ana *r q cosh^Z + a sinhqtfj w 


The roots inp are real and negative, and we can proceed to an inter- 
pretation by the partial-fraction rule as usual. Or, using the expansion 
in * waves/ we have 


If the length is great enough to make the terms involving e~ nl inap- 
preciable, we can reduce this to its first two terms, thus : 



If a is great, the solutions reduce to those of 5.2 ; this is to be expected, 
for (1) then implies that V « 0 when x = 0, which is the boundary con- 
dition adopted in 5.2. If a is small, V reduces to S; the reason is that 
this implies that there is no loss of heat from the end, and therefore the 
temperature does not change anywhere. For intermediate values we may 
proceed as follows. If 


Bromwich’s rule gives 



■(f) 


Put 


If a h A }$x\dX 


X = K 2 . 


( 8 ) 

O) 


The path of integration for k is a curve going from Re~^ wl to Be^ 1 , 
where R is great, passing on the way within a finite distance of the 
origin on the positive side. Denote this path by M. Then 


I dK 


exp ?) dK - 


( 10 ) 

( 11 ) 
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= 0-w 

(l2) 

The second part of (11) can be written as an integral with regard to /*, 
where 

I* - k + ah (13) 

along a path obtained by displacing M through a distance ah ; but the 
integrand is regular between these paths and the route M may still be 
used. The second part is therefore 

- ~ - exp | pH - /* (2akt 4* j exp ( a2 hH + ax) dfx. 


= - exp (aVbH + ax) 1 1 - erf } ( 14 0 

and ?y = l-erf f-exp^ + a#) {1 -erf (£ + y)}, (15) 

where f = x/2h$ ; y = ah$ (16) 

Hence V = $[e rf f + exp (y 2 + ax) {l - erf (£ + y)}] (17) 

This is the same as Biemann’s solution*. 

The temperature at the end is 

F^o~£expy 2 {l-erfy}, (18) 


whence the temperature gradient at the end follows by (l). For small 
values of the time the temperature at the free end falls continuously ; 
the temperature gradient there is not instantaneously infinite as in 
5.4 (5). For great values of the time we can use the approximation 2.4 
(28), giving 

F “ =o= a/770rf) 

This is equivalent to one found by Heavisidet. Heaviside gives also an 
expansion in a convergent series, suitable for small values of t, but it 
is probably less convenient than the equivalent expression (18) in finite 
terms. We see from (19) that the longer the time taken the better is 
the approximation to (3 V/dx)x-o given by the simple theory of 5.2, 

* Biemann-Weber, JPartielle Elfferentialgleicliungen , 2, 1912, 95-98. 
t Electromagnetic Theory , 2, 16. 
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5.6- A long rod is fastened at the end x = 0, the other end not being in 
contact with a conductor. Initially it is at temperature 0, but at time 0 
the clamped end is raised to temperature 8 and kept there. Each part of 
the rod loses heat by radiation at a rate proportional to its temperature. 

The differential equation is now 


a V a a v 



(i) 

where a is a constant. Let us put 


j P + a 2 = /iV 2 

(2) 

and write the equation 

3 a F OTr 

(3) 

-r 2 F-0 

dx 2 


At x = l there is no conduction out of the rod, and therefore 3 V/dx 
vanishes. Also F= 8 when x ~ 0. The solution is 


«> 

This can be expanded in powers of er 2rl . The first term, which is the 
only one we require, is 

V=Se-™ ( 5 ) 

8 ( ( (\ + a a )V| d\ ... 

h-fT <•> 

put ' A. + a a = k 2 (7) 

Then 

But if K = a + /4. ; ( 9 ) 

the term in 1 /(« - a) becomes 

<*« 

as in 5.5 (14). The complete solution is 
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When a = 0 this reduces to , 

K - s (>- erf Si) <**> 

so that the disturbance of temperature spreads along the rod, the time 
needed to produce a given rise of temperature at distance x from the 

end being proportional to x\ If t is small enough to make small, the 
error functions will be practically unity except where x is not great 

compared with 2 h$. For such values of x the exponentials are nearly 
unity. Thus at first the heating proceeds almost as in the absence of 
radiation from the sides of the rod. 

But if is great and x/2 h$ small or moderate, the first error func- 
tion in (11) is practically - 1, and the second + 1. Thus in these con- 
ditions 

F (13) 

This is seen by a return to the original equation to be the solution 
corresponding to a steady state. This will hold so long as a$ - x/2h$ 

is large and positive, even if x\2h$ is itself large. The region where the 
approximation (13) is valid therefore spreads along the rod with velocity 

2 ah ; if at- has once become large, each point on the rod reaches a nearly 
steady temperature at a time rather greater than x/2ah. 

If atfi is large and x/2 h$ still greater, the solution is nearly 

{ x - 2a hu 

i+erf ( 14 ) 


If further the argument of the error function is large, we can write 


1 

2 \Ar 


V = t ; r Sv~ aCClh exp 


{■ 


(,r- 2a Mf \ 2&* 

4 AH J x — 2 a lit 


w'-lh, <“> 


In the regions tliat have not reached their steady state, the temperature 
resembles that for the problem without radiation, except that the small 

factor must be introduced. 

" x - 2a ht 


5.7. The cooling of the earth . Cooling in the earth since it first became 
solid has probably not had time to become appreciable except within a 
layer whose thickness is small compared with the radius. It is therefore 
legitimate to neglect the effects of curvature and treat the problem as 
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one-dimensional. Radiation from the outer surface must have soon 
reduced the temperature there to that maintained by solar radiation, so 
that we may suppose the surface temperature to be constant and adopt 
it for our zero of temperature. The chief difference from the problem of 
5.2 is that we must allow for the heating effect of radioactivity in the 
outer layers. Suppose first that the quantity P defined in 5.1 is equal 
to a constant A down to a depth II, and zero below that depth. Take 
the initial temperature to be S+ mas, where m is a constant. Then the 
subsidiary equation is 

0 <m<E 

and + mas) II <as, 

and the solutions are 
A 

V = -Trys + jS? + mx + Be~u x + Ce<s x 0<as<H\ 

f (2) 

and = JS + mas + De~Q x II<as , J 

A term in e qso is not required in the solution for great depths, because 
it would imply that the temperature there suddenly dropped a finite 
amount at all depths, however great, in consequence of a disturbance 
near the surface. The conditions to determine B , G, and B are that V 
vanishes at # = and that V and dV/dco are continuous at as = 11. 


Hence 

B + G + S+A/p^O, (3) 

i?0-« 2Z ’ + Q e qH +A/p = Be~ qH , (4) 

Be~ ** - Cetx = Be- (5) 

Solving and substituting in (2), we find 

A A 

V= S (1 — 0“®*) + mas + ~ {1 — — 0-2# si n h ^j. o <qs<H 

A f («) 

V = JS(1 - + mas + — (cosh qH— 1) as>H. 

P 


These solutions involve terms of the form ~ &- qa , where a is positive. To 
interpret, we write 

-LL ("■'-?)* 
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*P 




on integrating by parts. But 

= ( 8 ) 

and a further integration by parts gives 

L t) * ■ - hi (” ■ - e) ' “P ( rt - t) * 

...(9) 

InaH |r~- (.+£) (l - erf^j)- f v'j'--"- "■(“> 

Thus we can obtain a solution in finite terms, which is easily seen to be 
identical with that previously obtained by a more laborious method *. 

The temperature gradient at the surface is found most easily from 
(6) ; we have 

(ll) 

Bat i*-«-*s/ i eip(M-*Tr) J < 12 > 

h [ f 2 , kH\ d>K 

r'"r) 7 

= 2hHqe-« B ~He-i H 

= 2 h (tjirf r *W» * _ H (l - erf J~) (13) 

and therefore 

(14) 

This is identical with the result given in earlier works of mine except 
that in these the factor 2 was omitted from the last term, which is small. 

In the actual problem a considerable simplification arises from the 
fact that II is small compared with 2h$. On this account we can 
expand the solutions in powers of H and retain only the earlier terms. 
Thus for the surface temperature gradient we have 

CsrL 

. « + « + ^(i- * ■ ) (15) 

k (jt^ ^ ' 2k (yty 

* Jeffreys, Phil. Mag. 32, 1916, 676-691; The Earth, 1924, 84. 
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and for the temperature at depths greater than H 

■ V= $(1 - e' iX ) + mx + . e~ ax 

, /„ AH\ „ x AH* 

=’“ + ( 4, -ir; erf si + iF <“) 

5.71. An alternative possibility is that the radioactive generation of 
heat, instead of being confined entirely to a uniform surface layer, may 
decrease exponentially with depth. In this case the subsidiary equation is 

1 ¥-q‘V=-j a e-’ x -q t (S+mx) (1) 

at all depths. We already know the part of the solution contributed 
by S + mx. The remainder is 

r fr _Ae- aa -e-^ 

W h* q*— a* ( 2 ) 

q—om g-ax 1 

Bn * ( 8 ) 

and - = — 

q 2 -a? * \q + a q-aj 

= \ [e -9 ® - exp (y 2 + ax) {1 - erf (£ + y)}] 

+ i [e -°* -exp (f-ax) {1 — erf ($- y)}], ( 4 ) 

w ^ lere £=xl 2 h$ ; y = ah$. (5) 

Hence JF=-^- 2 {exp (y 2 - ax) - exp (- ax)} 

A 

+ £v [(1 - erf £) - £ exp (-f + ax) {1- erf(| + y)} 

A - 1 exp if -ax) {1 - erf(l-y)}] 

= Wa* [l-erf^-r‘ l "*-^exp (y 2 + ax) {1 - erf (£ + y)} 

+ iexp (f-ax) {1 + erf(y-£)}], ...(g) 
which is the same as the solution obtained by Ingersoll and Zobel*. 

The contribution of radioactivity to the temperature gradient at the 
surface is 

/dW\ _ A 1 

\ 9 a ) Xss0 Wq + a (0 

= Fa {l ~ eXp ^( 1-erf V)} (8) 

Mathematical Theory of Heat Conduction , Ginn, 1913. 
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When y is great, as it actually is, we have 

(Hr ^( 1 ~ yx) = (*’“ «* -s/(*o) ' 



5.8. The justification of the method is easier in problems of heat 
conduction than in those of the last chapter, because the integrands 

always contain a factor exp ^ — ^). This tends to zero when A, tends* 

to c±l co in such a way that the integrals obtained by differentiating 
under the integral sign always converge, and can therefore be sub- 
stituted in the differential equation directly. But the integrals for the 
temperature are of the form 

r.i//w«e(«-4*)£ 

and when we substitute in the equation 


dV 

dti 


a a F 


the integrand vanishes identically. 


CHAPTER VI 

PROBLEMS WITH SPHERICAL OR CYLINDRICAL 
SYMMETRY 

6.1. So far we have treated only problems of wave transmission 
or conduction of heat in one dimension. If our system has spherical 
symmetry, the equation of transmission of sound takes the form 

~~c 2 V^ = 0, (1) 

where $ is the velocity potential, and 

1 3 /_.»\ 


2 3$ 
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The differential equation is therefore equivalent to 


0 s 

c‘df 


(r*)-p(r*)-0. 


( 3 ) 


This is of the same form as the equation of transmission of sound in 
one dimension, r$ taking the place of but differences of treatment 
arise from differences in the boundary conditions. A similar transforma- 
tion can of course be applied to the equation of heat conduction. 

If the sjnnmetry is cylindrical, we take Cartesian coordinates y } z, 
and put 

a? + y* = w 2 , y — cc tan <j>, (4) 

and if & is independent of z and </> we have 


V 2 4> = 


ijw ?*\ 

vr dvr \ dvr) * 


( 6 ) 


which is capable of no simple transformation analogous to that just 
given for the case of spherical symmetry. This fact is connected with 
striking differences between the phenomena of wave motion in two and 
and three dimensions. 


6.2. Consider first a spherical region of high pressure, surrounded 
by an infinitely extended region of uniform pressure; the boundary 
between them is solid, and the whole is at rest. Suddenly the boundary 
is annihilated. Find the subsequent motion. The problem is that of an 
explosion wave. We suppose the motion small enough to permit the 
neglect of squares of the displacements. At all points 

>• <*) 

Initially there is no motion, so that is constant, and may be taken 
as zero. The excess pressure P is ~pd$/dt f where p is the density; 
this is initially a positive constant P 0 when r < a , but zero when r > a. 
Then we can take the subsidiary equations to be 

• = ? ' <a l (2) 

= 0 r>a.J 

The pressure must remain finite at the centre. Hence 

+ a sinh^ r<a (3) 

pp o v ' 

= B exp (- prjc) r> a (4) 

There can be no term in exp (pr/c) in (4) since it would correspond to 
a wave travelling inwards. Now the pressure and the radial velocity 


i 


ft 
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must be continuous at r = a ; hence $ and d&/dr must be continuous. 
These give 

( 5 ) 

- — 4- - J. cosh^ = --B exp (-— ) , (6) 

whence 

- 4 =^(c + op) B=~±(c-ap) eP *!° (c + op) e-J* 1 / 0 . 

(7) 

Thus outside the original sphere 

T ^ = ^ " «P) (r " a)/c “ ^2 + ap) e-*<r+*V*. . . .(8) 

The associated pressure change is 


p= 

-if 175 

2 r L\J0 

■ - aj e~P (•■-»)/#— Q + e~P . 

...(9) 

But 

c 

— a 
P 

= 0 when t is negative [ 

= ct -a when t is positive J 

....(10) 

Hence 

and 

fc \ 

— a 
\p > 

) g-p {»■-«)/« = 0 when ct < r - a 
-ct~{r- a)- a 
= ct-r when ct> r — a. 

....(11) 

Similarly 

d + “) 

1 e -v (r+a)/o _ q w ]^ en cif < r + a 


and 

= ct~r when ct> r + &. ... 

....(12) 


Hence the pressure disturbance is zero up to time (r - a)/c, when the 
first wave from the compressed region arrives, and after (r + a)/c, when 


the last wave passes. At intermediate times it is equal to ~ (r 


■ Gt ). 


Thus it is equal to P^aj2r when the first wave comes, — P Q a/2r when 
the last leaves, and varies linearly with the time in between. The com- 
pression in front of the shock is associated with an equal rarefaction in 
the rear*. 

Within the sphere the pressure is 

^ — [*- - ~ a (c + ap) sinh — 0 -^/°"] 
r Lp p 1 v c J 

= p t 1 ~ h- (| + a ) ( e ~ p {a ~ r)l ° ~ e ~ pia+r)l °)\ • • • -(i3) 


* Of. Stolces, Phil. Mag. 34, 1849, 62. 
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This is equal to P 0 up to time ( a-r)/c, then drops suddenly to 
P 0 (1 - a/2r), decreases linearly with the time till it reaches - Phifer at 
time ( a + r)/c , and then rises suddenly to zero. The infinity in the 
pressure at the centre is only instantaneous, for the time the disturbance 
lasts at a given place is 2 rjc, which vanishes at the centre. It is due to 
the simultaneous arrival of the waves from all points on the surface of 
the sphere; at other points the waves from different parts of the surface 
arrive at different times, giving a finite disturbance of pressure over a 
finite interval. If r<\a , the pressure becomes negative immediately on 
the arrival of the disturbance. 

The behaviour of the velocity at distant points is similar to that of 
the pressure. If u is the radial velocity, 

u - d$>jdr 

= -p$/c-<1>/r (14) 

If r is great the first term is proportional to 1/r, the second to l/r\ 
The first is therefore the more important. But the first term is simply 
a multiple of the pressure. Hence it gives no motion at a point until 
time (r~&)/c, when the matter suddenly begins to move out with 
‘velocity P 0 a/%pcr. This velocity decreases linearly until time (r -b a)/c t 
when it is - Poafepcr, and then suddenly ceases. The total outward dis- 
placement contributed is zero. The second term, however, gives a small 
velocity which vanishes at the beginning and end of the shock, and 
reaches a positive maximum at time rjc. It produces a residual displace- 
ment, of order ajr times the greatest given by the first term; this 
represents the fact that the matter originally compressed expands till 
it reaches normal pressure, and the surrounding matter moves outwards 
to make room for it. 


6.21. Consider next the analogous problem with cylindrical sym- 
metry. With analogous initial conditions, the subsidiary equation is 


-zar dm \ dm/ c 2 




The solutions are Bessel functions of imaginary argument and order 
zero, I Q (pm/c) and Ko^pm/c). The latter is inadmissible within the 
cylinder, because it is infinite when -sr^O. The former cannot occur 
outside it, for the following reason. The interpretation is to be an 
integral along a route through values of the variable with positive real 
parts, and when rs is great the asymptotic expansion of I Q contains 

exp (A m/c) as a factor. Hence the solution would involve exp (par/c) 
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and therefore a disturbance travelling inwards. The solutions are therefore 

* = -— 1 +AlJZZ) 0 <®<« 

PP \c ) 


Also d<3?/dt and d®Jdm must be continuous at isr= a. Hence 


••• 

( 3 ) 

Ai;(ef) -J9X.'(«). ... 

00 

Also, with the definitions below, we have the identity 


(z) ic, (z) - /„ (s) (z) = 2/7 rz. ... 

( 5 ) 

Hence we find for points outside w = a,, 




. Is- tt) T 

(6) 

But /„ (z) = — j* exp (z cos 0)d6, 

( 7 ) 

2 r 

K {) ( z ) '= - / exp (- z cosh v) dv, ... 

7 r Jo 

(8) 

whence 


h = ~wJJ 0 fo x ex P x (' + « U08<J c cosh v ) 

cos 0 dkdtidv. 


O) 


Performing first the integration with regard to A, we obtain a function 
of the form II {t — b). Thus 

y ~ ™ l f cos 6d6d% (10) 

I\ ttcJq Jo k j 


where the range of integration is restricted by the condition that 

ct + a cos 6 - m cosh v > 0 (11) 

It follows at once that there is no movement at a place until time 
(•nr -d)lc. Integrating next with regard to v, we see that the admissible 
values range from 0 to cosir 1 {{at + a cos 0)jw} 7 provided the quantity 
in the parentheses is greater than unity. Hence 

e —a. ( m ,^ (*£!«»•)<». 

r, no J \ vs / 



( 12 ) 
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So long as ct>vr-a, there will be a finite range of integration with 
regard to 0; the lower limit is then always zero. If ct >'& + <%, the 
inequality is satisfied for all values of 0 up to 7r, and therefore ?r is the 
upper limit. If ct < w + a, the inequality is not satisfied when 9 = ?r, and 
the upper limit is cos'" 1 (m — ct)!a. The disturbance at any point may 
therefore be divided into three stages, the first until ct = m — a, the 
second from then till ct = m + a, and the third later. 

We are concerned chiefly with the pressure. This remains constant 
until ct = cr - a\ in the next stage 

P _ a rcos-*(vr-ct)la COS 6 do ( v 

Po v Jo {(ct + a cos 9) 2 ~w 3 }& 

and in the last a similar integral with the upper limit replaced by tt. 
By applying the transformation 

ct + a - us = 26 ; ct + a cos 9 - -sr = 2# cos 2 \j/ (14) 

and integrating on the supposition that b is small, we find that soon 
after the arrival of the wave 

W©(H;-ia < I5 > 

When the wave arrives the pressure therefore jumps to W'®') , and 

then falls by-^+ of itself per unit time. The corresponding fraction 

in the spherical problem is c/a. At time -zv/c the pressure is still positive ; 
but when ct = -zzr + a or more and 9 is greater than the integrand in 
(13) is numerically greater than for the supplementary value of 9 , and 
thus P is negative. It actually tends to - <x> at ct = -a + a, and returns to 
finite negative values for greater values of t. The passage of the wave 
of rarefaction is therefore indefinitely protracted. To find out how it 
dies down with the time let us suppose that ct is greater than w, and 
that a is small compared with either. Then 

a?ct 

P = (c 2 1 2 - -nr 3 / 

approximately. The residual disturbance falls off like t~ 2 . 

6.22. If the motion was one-dimensional, as for instance if the original 
excess of pressure was confined to a length 2a of a tube, the resulting 
disturbance of pressure would consist of two waves, each with an excess 
of pressure equal to travelling out in opposite directions with 
velocity c. Comparing the results for the three cases, we see that the 
first disturbance at a given point outside the region originally dis- 
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turbed, in each case at the same distance from the nearest point of it, 
occurs in each case at the same moment. The increase of pressure in the 

one-dimensional problem is |P 0 , in the two-dimensional one |P 0 (a/w)^, 
and in the three-dimensional one $P 0 (a/r). In the first case the pressure 
remains constant for time 2a/c, and then drops to zero and remains 
there. In the cylindrical problem it begins to fall instantly, and becomes 
negative in an interval less than 2 a/c] it tends to - go , and dies down 
again asymptotically to zero. In the spherical one it decreases linearly 
with the time and reaches a negative value equal to the original positive 
one at time 2 a/c ; then it suddenly becomes zero again. 


6.3. Diverging waves produced by a sphere oscillating radially *. 
Suppose that a sphere of radius a begins at time 0 to oscillate radially 
in period 2? r/n. We require the motion of the air outside it. 

The velocity potential 3> satisfies the equation 

CD 

-Initially all is at rest ; the solution is therefore 

r®=zAe~ vrl \ ( 2 ) 

When r~a the outward displacement is, say, - sin nt when t> 0, and the 

n 

outward velocity cos nt. Hence 


and 


?•<£ = — 


EML' ® 

(f‘ j + w a ) (I +pa/c) 6X ^ { c^ 7 a ^\ 


. 3^3 


<ra 


6 ,a + d l n l 

;.2 


exp 


H (r - s) } {• 


t na . 
cos nt + ~ sin nt - exp 


(-!)} 


<>a 3 r f r-a } na . f. r-a] 

+ ® J » J L l c ) c l c ) 


•w 


when ct>r~a. 

The solution has a periodic part with a period equal to that of the 
given disturbance, together with a part dying down with the time at a 
rate independent of ?r, blit involving the size of the sphere. As there 
is no corresponding term in the problem of 6.2 we may regard it as 
a result of the constraint introduced by the presence of the rigid sphere. 

* Love, Froc. Lond. Math . Soc, (2) 2, 1004, 88 ; Bromwich, i&. (2) 15, 1916, 481. 

J 6 


* 



i 

I 
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Its effect on the velocity or the pressure is to that of the second term 
in a ratio comparable with {cjndf. 

B .4 ■ A spherical thermom eter bulb is initially at a uniform temperature 
equal to that of its surroundings. The temperature of the air decreases 
with height, and the thermometer is carried upwards at such a rate that 
the temperature at the outside of the glass varies linearly with the time. 
Find how the mean temperature of the mercury varies*. 

The temperature within the bulb satisfies the equation 


w(rV)-q'rV= 0, (1) 

■where p = Wq* (2) 

That at the outer surface of the glass is Gt, where G is a constant. 
But the glass has only a finite conductivity, so that the surface condition 
at the outside of the mercury is 

. d l=K(Gt-V).. ( 3 ) 

K being another constant. The solution of (1) is 

F=~sinh qr (4) 

and (3) gives A = ^ m ^ KG/p — (5) 

w Ka smh qa + qa cosh qa - smh qa ’ v J 


where a is the^inner radius of the glass. 

The mean temperature within the bulb is 


i=- a /v i 

Jo 


= (qa> cosh qa - sinh qa) 

_ 3 KG qa cosh qa - sinh qa 

apq 2 Ka sinh qa + qa cosh qa - sinh qa 

In applying the partial-fraction rule, we notice that near p = 0 

v = 3 - gg (fr + m 

0 qa(Ka + -J* KfcP + J (fa?) 




Bromwich, Phil, Mag. 37, 1919, 407-41 9 j A. B. McLeod, Phil. Mag. 37, 1919, 134. 



PROBLEMS WITH SPHERICAL OR CYLINDRICAL SYMMETRY 83 


4V 




V' 




so that there is a constant lag in the temperature of the mercury in 
comparison with that of the air. 

The other zeros of the denominator give exponential contributions, 
which are evaluated in Bromwich's paper. 


6.5, A cylinder of internal radius a can rotate freely about its axis. 
It is filled with viscous liquid, and the whole is rotating as if solid with 
angular velocity o> 0 . The cylinder is instantaneously brought to rest at 
time t-0 i and immediately released. Find the angular velocity later* 
(Math. Trip. Schedule B, 1926.) 

The motion is two-dimensional, and there is a stream-function ^ 
satisfying the equation 


( V V-’) V>. 


9 (».y) ’ 


(i) 


where v is the kinematic viscosity. Since the motion is symmetrical 
about an axis the right side is identically zero. Put 


l=p, p-n>, V’=i^( w ±.) (2) 

Initially \j/ = £ vr a) 0 ; V 2 ^ = 2w 0 * (3) 

and the subsidiary equation is 

(V 2 — r 9 ) V 2 ^ = - 2 r 2 o) 0 (4) 

The solution is 


= AIq {rm) + BKq (?'*zzr) + G log or + D 4- <d 0 (5) 

The velocity must be finite on the axis ; hence B and G are zero. Also 
D is independent of ra* and therefore cannot affect the motion. 

If / be the moment of inertia of the cylinder per unit length and <*> 
its angular velocity, the equation of motion of the cylinder is 


= -Snaps, ( 6 ) 

where p H is the shearing stress in the fluid. Now 

/3V av\ 

=vpa L(l*£) w 

evaluated on the outer boundary at the point ( a 9 0). Hence 

p, = vp~r [ral 0 (pa) - 2 1 0 ' (ra)J (8) 

Cl 


6-2 
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Since the cylinder starts from rest, 

ipoj = - 2 ? rvp Ar [ra I 0 (ra) -2 /</ (ra)] (9) 

Also a<o must be equal to the velocity of the fluid. Hence 


aw =Ar /<>' (ra) + aw Q (10) 

Eliminating A we have 

a> \_{Kp - 2) / 0 ' (ra) + ra / 0 (ra)] = <*> 0 [ra / 0 (ra) - 2/ 0 ' (ra)], ...(11) 

where %TrvpaK=I, (12) 

The operational solution is therefore 

mJ 0 (?-a)-2/ 0 '(m) , 

0 (Kp - 2) 4' (m) + ra /, (ra) ( j 

But 

/ 0 (»*«) = 1 + ir'a 2 + + ... ; /<,' (ra) = £ra + T V^« 3 + .... (14) 

The contribution from =0 is found to be 

<o 0 a a „ ^ 

a- + 4vjST ~ < * >1 ’ ^ 

say. This is the ultimate angular velocity. The other terms arise from 
the zeros of the denominator. If we write Jc for r, and substitute for K 
in terms of ah,' we find that these satisfy 


|(o> 0 - wj) + 2 j J 1 (Jed) - JcaJ n (Jca) = 0, 

The solution is then of the form w x + The coefficients A k are 

determinable by the usual method. 


CHAPTER VII 


DISPERSION 


In the propagation of sound waves in air and of waves on strings the 
velocity of travel of waves is independent of the period. In many problems 
this is not the case ; waves on water afford an important example. 


7.1. Consider a layer of incompressible fluid of density p and depth 
ZT. Take the origin in the undisturbed position of the free surface, the 
axis of z upwards, and those of w and y in the horizontal plane. Put 


d 

a r* m > 




a 2 a 2 a 

— j- — — a.- 

da? dy* 


■CD 
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The velocity potential satisfies the equation 

9 2 3 > 

— + = 0 ( 2 ) 

At the bottom the vertical velocity vanishes. At the free surface, so long 
as the motion is only slightly disturbed from rest, we have 


3s ~ dt ’ 


( 3 ) 


where £ is the elevation of the free surface. Then we must have 
s, p oosr(s+jr) 
r sin rH 




The pressure just under the free surface is -Tp?^, where Tp is the 
surface tension. But by Bernoulli's equation it is also equal to 

-pp$-gpt+F(t\ 

and F (t) does not affect the motion and can therefore be omitted. We 


have therefore the further surface condition 

- p$ = (g- Tr 2 )& (5) 

Combining this with (4) we have the differential equation for £ 

{p 1 ~(g~ Tr*) r tan rH) £ = 0 (6) 


If the fluid starts from rest with £ equal to to, a known function of x and 
y i the subsidiary equation will have a term on the right, and the 
operational solution will be 

^ = j» a — (i? — TV) r tan rH 

= cosh t {( g - TV' 2 ) r tan rH$. £ 0 (8) 

In the corresponding problem of a uniform string in 4.1 the coefficient 
of t was simply pc. This solution, like d’Alembert's solution for the string, 
assumes that there are no ends to reflect waves. 

Suppose first that the original disturbance consists of an infinite 
elevation along the axis of y, with no disturbance of the surface 
anywhere else. Then £ 0 can be replaced by DH{x) i and (8) is equivalent 
to the integral 

£ = [ eft* cosh t {(g — TV) k tan kH}% d k (9) 

27TI J Jj 

The integrand has an essential singularity wherever kH is an odd 
multiple of L cannot therefore cross the real axis within a finite 
distance of the origin, but becomes two branches extending to + oo above 
and below the axis. 
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7.2. The method of steepest descents. The elevation of the surface is 
thus expressed in terms of integrals of the type 

S= exp t{f(e)}dz (1) 

where f(z) is an analytic function, and t is real, positive, and independent 
of z. Put, following Debye, 

f(z) = R +//, (2) 

thus expressing it in real and imaginary parts. If the integral is taken 
along an arbitrary path, the integrand will be the product of a variable 
positive factor with one whose absolute value is unity, but which varies 
in argument more and more rapidly the greater t is. There will evidently 
be advantages in choosing the path in such a way that the large values 
of R are concentrated in the shortest possible interval on it. Now if 


z^a + iy, ( 3 ) 

we shall have _ + _ = 0 ; -^ = 0 (4) 

It follows that R can never be an absolute maximum. But it can have 
stationary points, where 

dR dR 

(5) 


and we know that these points will also be stationary points of /and 
zeros of /' (z). These points are usually called the * saddle-points/ or 
sometimes * cols.' Through any saddle-point it will in general be possible 
to draw two (sometimes more) curves such that R is constant along 
them. In sectors between these curves R will be alternately greater and 
less than at the saddle-point itself. The sectors where R is greater may 
be called the 'hills/ those where it is less the 'valleys/ If our path of 
integration is to be chosen so as to avoid large values of i?, it must 
avoid the hills, and keep as far as possible to the valleys. If then the 
complex plane is marked out by the lines of R constant through all the 
saddle-points, and A and B lie within the same valley, our path must 
never go outside this valley ; but if A and B lie in different valleys, 
the passage from one valley to another must take place through a saddle- 
point. In the latter case the value of the integral will be much greater 
than in the former, and therefore interest attaches chiefly to the case 
where the limits of the integral lie in different valleys. 

The paths actually chosen are specified rather more narrowly; the 
direction of the path at any point is chosen so that | dR/ds | is as great 
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as possible. If ^ is the inclination of the tangent to the path to the axis 
of x, we have 


9 R ,m . 3 j? 

^ = cos^- + sin^ ¥ , 


( 6 ) 


and if this is to be a numerical maximum for variations in if/ 


a • . dR 
0 = - sm t/f ^ + cos 


dR 

dy 


dB = dJ 
dn 9 s 9 


00 


where dn is an element of length normal to the path, drawn so that 
when ds is in the direction of x increasing, dn is in the direction of y 
increasing. Hence I is constant along the path. Such a path is called 
a ‘line of steepest descent/ There will be one in each valley. In general 
the limits of the integral will not themselves lie on lines of steepest 
descent, but can be joined to them by paths within the valleys. 

In general lines of R constant through different saddle-points will 
not intersect ; and there will be only one saddle-point on each line of 
steepest descent. For the former event would imply that R has the same 
value at two saddle-points, the latter that I has, and either of these 
events will be exceptional. It follows that as we proceed along a line of 
steepest descent R will rarely reach a minimum and then proceed to 
increase again. For if R had a minimum 9 R/ds would be zero • but 9 1/ds 
is zero by construction, and therefore the point would be another saddle- 
point. Lines of steepest descent usually terminate only at singularities 
off(z) or at infinity. 

The path of integration once chosen, the greater t is the more closely 
the higher values of the integrand will be concentrated about the saddle- 
points. Thus we can obtain an approximation, which will be better the 
larger t is, by considering only the parts of the path in these regions. 
In these conditions we can take 


/GO =/(*<>) + i (* “ *<>)*/" (*»)> (8) 

where z 0 is a saddle-point. Put 

l/Whi] |*-*|=»r. (9) 

Then on a line of steepest descent 

/(*) =/(*o) - idr 2 , 


( 10 ) 
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on the side after passing through z 0 . Then 

B-*f exp {tf(z Q )} exp (- ^Atr 2 ) d (r exp to) 

= P {</(*») + «*}i (12) 

the approach to z 0 and the descent from it making equal contributions. 

To get from A to B it may be necessary to pass through two or more 
saddle-points, with probably traverses in the valleys between the lines 
of steepest descent; Then each saddle-point will make its contribution 
to the integral. 

The error involved in this approximation arises from the terms of the 
third and higher orders omitted from (8). Its accuracy therefore depends 
on exp (— having become small before exp (s 0 )|) has 

begun to differ appreciably from unity. Hence 


(13) 


must be large. In most cases the approximation is asymptotic, and does 
not represent the first term of a convergent series. 


7.3. In problems of wave motion we often have to evaluate integrals 
of the form 

. £ = hi Izf W exp («*■ -• yt)dx, (1) 

where (*) and y are known functions of *. As a rule ( k ) is an even 
function, and y an odd one. When k is purely imaginary y is also purely 
imaginary. We require the motion for large values of t, and possibly 
also of x. 

The function 4> (/c) usually introduces no difficulty. It does not involve 
x or t, and therefore when these are large enough it can be treated as 
constant throughout the range where the integrand is appreciable. 

It is usually convenient to replace * by t k and y by tyj and to consider 


the equivalent integral of the form 

£ = J ^00 exp i(i«0-yO<fc, (2) 

The saddle-points are given by 

( 8 ) 


so that a given ratio xjt specifies a set of predominant values of k 0 . But 
y' is an even function of k , and therefore if k 0 is a saddle-point, - k 0 will 
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be another, and if the adopted path passes through either it will pass 


through the other. We may take k q real and positive. Also 

/" (*o) =- iyo"*< (4) 

Thus if y 0 " is positive k- k q has argument - J tt on the line of steepe 
descent, and the contribution from k q is 

5^ (j^{) ^ (*•) e ~ iLW ex P 1 ( K ° x ~ r» # ) ( 5 ) 

The contribution from - k 0 (where y" is negative) is similarly 

V' («<D e i ‘ x exp {- 1 (K 0 a? - y 0 «)} , (6) 

and the two together give 



Similarly if y 0 " is negative the two saddle-points give 

/ 2 

'I' 0<o) COS (k 0 X -yot + lir) (8) 

7.31. These formulae, due to Kelvin, are the fundamental ones of the 
theory of dispersion. Consider first the cosine factor, and suppose x in- 
creased by 8x and t by 8t. Then ^x-y Q t is increased by 
k 0 8# - y 0 8£ + {x - y 0 r t) &c 0> 

the term in 8k 0 appearing because * 0 is defined as a function of x and t 
by (3). But the coefficient of 8 k 0 is zero by (3). If then t is kept constant, 
£ will vary with x with period 2 tt/k 0 ; and if x is kept constant, £ will 
vary with t with period 27r/y 0 . Hence 2 ? t/k 0 is the wave-length, and 
27 r/y 0 the period, of the waves passing a given place. A phase occurring 
at a given place and time is reproduced after an interval at a place 
such that 8x = y 0 8t/K Q . Hence y Jk 0 is the velocity of travel of individual 
waves. It may be denoted by c, and called the wave-velocity. 

But k q has been defined by the equation 

x-y Q '£ = <), 

so that a given wave-length and period always occur when x\t has a 
particular value ; they seem to travel out with velocity y 0 ', which is 
called the group-velocity. It may also be denoted by <7. In general the 
wave-velocity and the group-velocity are unequal, so that a given wave 
changes in period and length as it progresses. They are evidently con- 
nected by the relation 


( 9 ) 
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7.4. Keturning now to 7. 1 (9) we can separate the hyperbolic function 
into two exponentials, which will represent wave-systems travelling out 


in opposite, directions. One of them is equivalent to 

£ = r- J exp t (kos - yt) (Ik, (1) 

47T y -oa 

where y 2 = (g + 2V*) k tanh k i?. (2) 

When k is small, 

+ ( 3 ) 

+ W 

When « is great, 

c = (rK)i; C«$(ZV)* (5) 


In all ordinary cases Tjg is insignificant in comparison with IP. Hence 
for some intermediate value of * the group-velocity is a minimum; it 
tends to infinity for very short waves, and to a finite limit for very long 
ones. Three cases therefore arise. If acjt is less than the minimum group- 
velocity, there will be no saddle-point on the real axis, and the disturb- 
ance will be small*. If it lies between this minimum and (gtl)^, two 
(positive) values of k will give saddle-points, and each will contribute 

to the motion. If it is greater than {gUf 1 , the only saddle-point will 
correspond to a short wave. The disturbance at a given point will there- 
fore be in three stages. In the first, leading up to time xKglpft, only 
very short capillary waves will occur. Then long gravity waves will 
arrive, the wave-lengths of those reaching the point diminishing as time 
goes on. Superposed on them are further capillary waves, their length 
increasing with the time. At a certain moment the wave-lengths of the 
two sets become equal. This corresponds to the arrival of the waves 
with the minimum group-velocity. From then on the water is smooth. 

7.41, Two typical cases therefore arise according as the wave-length is 
large or small compared with the one that gives the minimum group- 
velocity. Take first gravity waves, such that aft is small compared with 


* It can be shown that the saddle-points are so plaoed that the relevant one 
contributes an exponential with a negative index to the solution. This is almost 
obvious from considerations of energy. 
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(gl-pft but large compared with the minimum group-velocity. In these 
conditions we can write simply 

f = g K ; c = (g/K)i; G^Ug/K)i; dC/d* — Kfflrf. ...(l) 
The solution is then 

(irf) 2 (^) « os ( ,c ° a, -y«* + i ,r ) C 2 ) 

wteM I ~ I (£>* < 8 > 

and the amplitude increases towards the rear of the wave train like 
aT* 

7.42. Take next the capillary waves, short enough for gravity to be 
neglected. Then 

y 2 = 7V; c = (2"k)*; 0=%(Tx)i; dC/d^UT/ic)* (1) 

At a given instant the amplitude is therefore proportional to or to 

The front of the disturbance is therefore composed of a series of 
capillary waves whose amplitude tends to infinity, and the time taken 
for them to arrive is infinitesimal. 

This impossible result arises from the form assumed for the original 
displacement. In taking f 0 ~DH(x) we assumed that unit volume of 
liquid was originally released on unit length of the actual line x = 0 in 
the surface. The mean height of this mass of liquid was therefore infinite, 
and its potential energy also infinite. The system being frictionless, this 
energy must be present somewhere in the waves existing at any instant, 
and infinite amplitudes are therefore a natural consequence of the initial 
conditions. If instead we suppose that the same volume of fluid was 
originally raised, but that it was distributed uniformly between x~± l, 
its elevation was 1/2/ in this range. Expressed in operational form this 
gives 

^ = ( 2 ) 

The appropriate solution can be found from 7.4 (1) by introducing a 

factor ( &lKl " e ~ lKl ) or — y ■ into the integrand. If the solution already 

found makes k 0 1 small, this additional factor will be practically unity, 
and the same solution will hold. Waves whose length is large compared 
with the extent of the original disturbance will therefore not be much 
affected by its finiteness. 
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But if k q I is large we must consider separately the contributions from 
the terms in e ld and e~ tKl . The factor 1/kI will give an extra 1/kJ in all 
the solutions. The result will be that at a given instant the amplitude 

of a gravity wave will vary like k 0 ~^ and that of a capillary wave like 
Waves whose length is short compared with that of the original 
disturbance are therefore heavily reduced in amplitude. 

On deep water the minimum group-velocity is 18 cm. /sec., corresponding 
to a wave-length of 4*6 cm. and a wave-velocity of 28 cm./sec. If the 
original disturbance has a horizontal extent of 1 cm. or so, only waves 
with lengths under 1 cm, or so will be affected, and the amplitudes of 
both gravity and capillary waves will increase steadily with diminishing 
group-velocity. A wave of large amplitude will therefore bring up the 
rear, and will leave smooth water behind it. This is observable in the 
waves caused by raindrops and other very concentrated disturbances. 
But if the extent of the original disturbance exceeds a few centimetres 
the capillary waves produced will be very small, and the largest 
amplitude will be associated with a wave whose length is comparable 
with the width of the disturbed region. The largest wave produced by 
the splash of a brick, for instance, has a length of the order of a 
foot. 

7.5. Two exceptional cases may arise in the treatment of dispersion, 
which are both illustrated in the present problem. The validity of the 
approximation 7.2 (12) depends on exp being proportional to 

exp At?*) on a line of steepest descent. If f"(z) has varied by a 
considerable fraction of A before this expQiiential has become small the 
approximation will not be good. This may happen if A is itself small, 
or if there are two saddle-points close together. Instances occur when 
there is a maximum or minimum group-velocity, or if the group-velocity 
tends to a finite limit when k becomes very small. In the former case a 
value of xjt a little greater than the minimum group- velocity will givo 
two slightly different finite values of k 0 . In the latter * 0 may bo small 
enough for the proximity of - * 0 to affect the contributions of both. 

7.51. Since y is an odd function of /c, we may suppose that when k is 
small 

y = c 0 K~c a K 8 + 0 (k 5 ) (1) 

and then 7.3 (2) is equivalent to 

$ M ^p i (kw - Co *t + CzkH) (Ik, 


( 2 ) 
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There are saddle-points where 



If, as in 7.3, if/ (*) is constant or tends to a limit different from zero when 
k tends to zero, this integral is nearly 
1 f 00 

£ = ~ Jo ^ (0) cos {/c (# - c 0 t) + c 2 kH} dK 


1 r°° 

~ ~ *^C0) oos(vP-mv)dv, (4) 

where m=(c„t- x )/(e»0* (5) 


The integral involved here is called an Airy integral*. It is finite and 
positive, but not stationary, when ffla'O; it has a maximum when 
m - 1*28, and oscillates for greater values of m with steadily decreasing 
amplitude. For negative values of m it tends asymptotically to zero. 



* Airy tabulates J cob ^ir (17 3 — 771*1?) dv in Canib. Phil, Trans. 8,1849,598. The 

graph, given here is adapted from Airy’s table. The integral can also be expressed in 
terms of BesBel functions of order J. See Watson, 188-190. 
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Consequently the disturbance we are considering produces an immediate 
rise of the level of the water at all distances greater than c Q t, though 
this rise is very small at great distances. The maximum rise of level is 
where x is rather less than c^t, so that it has not travelled out from the 
origin so far as a wave with the limiting velocity would. The maximum 
is followed by a series of waves of gradually diminishing length and 
amplitude, merging ultimately into gravity waves of the deep-water type. 

7.52. In the case where there is a minimum group-velocity for a finite 
wave-length, let us, with a somewhat different notation from that used 
so far, denote the minimum group-velocity by y 0 ' and the corresponding 


values of k and y by * 0 and y 0 . Put 

( 1 ) 

y = y 0 + yo'^i + 0 + }yo f/ Ki + (2) 

Then 


1 f 00 

£= — I i/r(K)expt(K 0 ar-y 0 #)expi(Ki^-y 0 , K 1 #-|-yo"'Kz 3 ^)^«i-"(3) 

A™ J~k q 

with an analogous contribution from the negative values of k • lienee 
with the same type of approximation as before we shall have 

2 f°° 

£ — - / ^ Oo) cos fox — yo cos x — y 0 t — ^yo £) d/cx 

irJo 

= - if/ ( k 0 ) cos (kqX - yot) (&yo"t)-^ f cos (v 3 - mv) dv , ...(4) 

7T JQ 

where m = — — — , (5) 

Qyrtf 

The solution is therefore the product of a cosine and an Airy integral, 
the interval between consecutive zeros being much greater for the 
former than for the latter. In the neighbourhood of a point that has 
travelled out with the minimum group-velocity the waves have the 
corresponding period and wave-length, but their amplitude falls off 
rapidly towards the rear. In front of this point the amplitude increases 
for a while and then oscillates. 

In both these exceptional cases we notice that the amplitude associated 
with the critical velocity falls off only like the inverse cube root of the 
time, whereas in the typical case it falls off like the inverse square 
root. Hence the further the disturbance progresses the more will the 
waves with the critical group-velocities predominate in relation to the 
others. In waves on water, however, these phenomena are often modified 
by the greater damping effect of viscosity on short waves. 


BESSEL EUEOTIOITS 


95 


CHAPTER VIII 

BESSEL FUNCTIONS 

8.1. Consider the operator p~ n exp (1/p), where n is a real number. 
By our rules we can write it as 

p~ n exp s ^7^= a “7-T 
r= o n (r) r =0 n (r) 

oo j.n+r 

r=o n (r) H(n + r) ^ 

= ■#»(*), (2) 

say, when t > 0* By differentiating the series in powers of t we have 

~E n (t)=M,^(t) (3) 

The function E n (t) satisfies the differential equation 
d*y n-ldy y 

df t dt~i~° 

The other solution of this equation may be taken as 

Fh ^ = Jo n>)n (»•-«) = Jo n^-«) ® 

unless n is an integer. 

In the form of a complex integral, 

B'V-bL r 0Xp ( X * + D^’ ® 

where M passes round the origin on the positive side and proceeds to 

- co at the ends. If t is positive this definition is equivalent to the 
series (1) without restriction on n. 

Now if we replace t by (is?) 2 , 

co a^VM-ar 

E n (i^ 3 ) = (i #)" 2 c 

i \i ) r=0 n(r)n(w+r) 

= (!*)•/.(*), (7) 

where I n (z) is now the modified Bessel function, satisfying the differ- 
ential equation 

*i{, d £)-V+^y-0. (8) 
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Prom (6) 

IM-(k*Y n ~ f M exv(±W + {)~, (9) 

which is equivalent to Schlafli’s form*. Put also 

= * (10) 

Then I n (*) ~ /^exp %z (k + \) ~ (11) 

if we take for z the positive value of (2t)K Now when A. is real and 
positive k is real and positive ; put now 

K k+ «) =/x (12) 

Then * = /x + - 1)& (13) 

the positive sign being taken for /x positive and greater than 1. Then 
when /x is + co , k is + cc , and as k describes a large circle the arguments 
of k and ju. are nearly equal. Thus the path M for k is equivalent to the 
path M for /x, provided only that we make it cross the real axis at a 
value of /x greater than + 1. 

Then /» (?) = — [ exp (us) r- (14) 

whence, in operational form, 

In (0 = J—~ T (l*>) 

This can be expanded in descending powers of p , and therefore corre- 
sponds to a series of ascending powers of t which will converge for all 
values of t, real or complex, and may be taken as the definition of the 
function. If n is fractional we remove ambiguity as to the argument by 
making I n (t) real and positive when t is real and positive. The coeffi- 
cients are determined by (7) ; thus 

CO /l /\«+2r 

/»(*) = 5 y-7^%7 \ (16) 

w r=oII (?■) n (n+r) 

The relation (3) is equivalent to 

~{t n In{t)} = t»I n ..,(t), (17) 

a familiar recurrence relation. 


* Watson, Bessel Functions , 176-6 and 181. 
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If in (16) we replace n by ~n, we obtain I- n (£), which is easily seen 
to be another solution of (8) ; but if n is a positive integer II (~w + r) 
is infinite for all values of r from 0 to n- 1, so that the first n terms of 
the series for I~ n (t) vanish, and the remainder are identical with I % ($)■ 


8,2. Now consider 


= f exp (fit) t T~, (I' 

taken along any path in the n plane. We have 

J,d fj, d(fi\ ^ , 2N if f jSy Aft 2 + tp-ft~M 2 J 


But 

?(/j: 2 -l) + tfx-n 2 ^ 


( 0 * 3 - 1 )* ***\ + ± 

i{#i+ Oi*-!)*}*^* J 


so that the integral is the difference between the values at the limits of 

pi 1 - !)** + » &llt (4) 

If then the path N is such that this quantity has the same value at 
both ends, </> (if) will be a solution of the differential equation. Now the 
integrand in c j> (i t ) has branch-points at ± 1, and also at infinity unless 
n is an integer. If then N is a closed contour and does not go round 
one or both of ± 1, (t) is identically zero. If t is positive and N begins 

and ends anywhere in the second and third quadrants respectively, and 
passes to the positive side of either branch-point, (4) vanishes at both 
limits on account of the factor el 1 *, and </> (t) satisfies the differential 
equation. 


8.21. Let N be the loop if already considered. Then </> (if) is I n (t). 
We may contract the loop vertically till it becomes the two sides of the 
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real axis from - oo to + 1, with small circular arcs about ± 1, as ii 
Fig. 8. On the lower side from - oo to - 1, we can put cosh v 
(/x 2 ~ 1)^ = e‘" t7r sinh v, and this part of the integral gives 

— /V w <r <cosh V e~ nv dv~-J~ [ e mir e~ tC0Bhv e-‘ nv dv. ...(5) 

2? nj w 2m Jo ' 

From - 1 to - oo , on the upper side of the axis, we take similarl] 

/x = e™ cosh % (p 2 ~ 1)~ = e lir sinh v, and the integral is 

^ f Q—n in - g-i cosh v g-nv 

2m Jq 

i r 

The two parts together give - — sin mr J e~ tco * hv e- nv dv. 

When we traverse the small semicircle below -1, /x and /x-1 keej 
the same argument, but that of jtx + 1 increases from -ir to 0. Hence 

from —1 to +1 we can put /x = e“ t7r cos 0, (/x 2 - 1)^ = 0 *"^* sin 0, with ( 
increasing from 0 to m The corresponding part of the integral is 

— r e n™ e -t™*o e -niOdQ = jL r 6 Um0 e niS dd t 
2 t rj 0 2-ir J 0 

On approaching + 1, /x is real and positive, while (/x a -l)^ still has 
argument - On traversing the circle about +1 we increase the 

argument of (/x a — 1)^ by m Hence on the path from + 1 to — 1 we put 
/x = cos0, (/tx 2 - 1)^ = e~ l7r sin# with 6 increasing from 0 to 7r, and the 

corresponding part of the integral is — f e tcm0 e~ nt9 dd, The two last 

Jq 

1 r 

parts together give - / cos nOe tcm6 dO, and in all 

vr Jo 

4 (t) m i jT OOS no e im « $ M - i sin m £ ...(6) 

8,22, By 8.2 we may alternately take N to be a loop about - 1 alone, 
as in Fig, 9. 
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t — % 

* m 




j* 


The contribution to the integral from the lower side of the path is as 
before. That from the upper side is different, for in describing the circle 

about - 1 , ijl keeps the same argument, while that of 1 )^ increases 
by 7 r, so that on the upper side we can put y. = e~ lir cosh % (ju . 2 - 1 )^ = sinh v, 


1 f°° 

and the integral is - — J e mrr e tcoahv e nv dv. Thus in this case 

6(t)=- —f e mrr e- tc03hv e- nv dv- — f e ntir e^ tcoabv e nv dv 
w 27Ttjo zttlJo 

1 f°° 

/ g nur cosh nve~ tC0Bhv dv (7) 

mjo 


Now if in ( 6 ) we replace n by -n we have 

/_„(£) = - f cosn0e tcosO d0 + -sin9i7r f e nv e" iC08]lv dv f ...(8) 
o v Jo 

2 r 

whence I~ n (f) - / u (t) = - sin mr j cosh m e^ tC08hv dv 


= sinwr. K n (t), (9) 

say. Then in (7) we can write 

<*>(*) =-J (10) 

and we may take 

K n (t) = - f cosh twer i0Mhv dv ( 11 ) 

V J 0 


as our second solution of the differential equation even in the cases 
where (t) coalesces with I n (t), Then ( 6 ) and ( 11 ) express our two 
solutions in real form*. 

If t-*- 0 , K n (i t ) evidently oo , and therefore has a singularity at the 
origin. I n (t) and K n (t) are both essentially positive. 


8.3. Asymptotic Approximations . If t is large, we can replace 
0 -tcoHhv j n ^ !) aiHiuh*^ w hidi is always less than e~ t e~^ vit . Hence 

the important values of this factor are concentrated within a range of v 
of order rt Within such a range cosh raw may be replaced by unity 
provided n l jt is small. Then to a first approximation 

= l f°er* e -i vH do = l0-*. 

TT Jo 7T V t 

(1) 

* This notation for IC n corresponds to Heaviside’s. Most writers on Bessel func- 
tions drop the factor 2/ir; Watson adopts this 
unfortunate. 


practice though he recognizes it a s 

library , 

/r 
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If we put 


*■«-</; 5*"( 1 + ? + ? + -) 0 > 


substitute in the differential equation, and equate coefficients of powers 
of t , we find the coefficients in succession. We get 



4^ 9 -l 3 
1! 8* 


(4w 2 - 1 2 ) (4ra 2 -3 2 ) 
+ 2 ! (8#) a 



The series is formally divergent for all values of t, but if t is large the 
terms decrease to a minimum, and it can be shown that the error in 
stopping at any term is less than the next term. 

Again, if t is large the first integral in (6) can be replaced, subject 
to n*jt being small, by 


/„(,). (4) 

By a similar process to that used for K n (t) we can obtain an asymptotic 
formula 


r /jA __ e l f, 4 if - 1 2 , (4rc 2 - l 2 ) (4 n Q - 3 2 ) 

nW sj(27rt) l 1 1 ! M + 2! (8tf) a 



8.31. If n*/t is not small, cosh nv in (11) may become great before 
exp (- 1 cosh v) has become small enough to swamp it. Then we can 
replace K n (t) nearly by 


K n (t) = ^ J exp (m - 1 cosh v) dv = ~ J exp f(v) dv, ...(1) 

(2) 


say, and the integrand has a maximum where 
/' (v) = ?i — t sinh v = 0. 

There 


sinh v = n/t ; f" (v) = - t cosh v ; f(y) = n sinh" 1 ? - (if 2 + w 2 ) K . ..(3) 

t 

Thus K n (t)=*Jl ...(4) 

To obtain an approximation to I n (t) we use 8.1 (14). We have 

l/W! (^ji « 

where /(/*.) =/*,£- w log {p + (/a 9 - 1)^} (6) 

There is a saddle-point where 






O a -l)4 


-0j 
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tlwe ,-(l. $*, f W - ^ - \ ($ * l)* (8) 

/O) = («* + O 4 - » log {(l + + || (9) 

Since /" (/*.) is positive the path of steepest descent is parallel to the 
imaginary axis; we find 

fF- ...(10) 

If t is kept constant and n is increased, the dominant factor in K % (t) 
is (2 n/t) n e~ n , which tends to infinity; the dominant factor in I n (t) is 
e” (2n/t)~ n , which tends to zero. 

8 . 32 . We note that, if y x and y a are any two solutions of the differ- 
ential equation, and 

& (i) 

% = y" Vi - y»" yi = - ^ (yi yr - y* yi) = - ^ (*), ( 2 ) 

so that $ (t) is of the form Aft, where A is a constant. If we take 
^i = 4(0j \ and use 8.3 (3) and (5), we find A by consider- 

ing large values of t\ then 

In 00 K n (0 - I n (0 Kn (f) = 2/irt ( 3 ) 

8.4. The ordinary Bessel functions may be treated similarly, starting 
from the operator 

P~ n exp (- 1/p ) = ( J tf J n (2* 4 ), (1) 

or j — - j— =J n (t) (2) 

But when (2) is converted into a complex integral the branch-points are 
at ±i, and the discussions of 8,21 and 8.22 need modification. This 
function satisfies the standard Bessel equation 

<■*-*» 

and is fully treated in the standard works of Watson and of Gray, 
Mathews, and MacBobert. 

8 . 5 . In - physical iDroblems Bessel functions usually arise with 
arguments pm/c. The foregoing work is still applicable because when 
an operational solution is interpreted as a complex integral p is replaced 
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by a number with a positive real part. Then the factor in the integrand 
is replaced by one that still tends exponentially to zero at //. = - oo, and 
the integral expressions survive without alteration. The importance of 


K n jo) is indicated by the factor exp 



in its asymptotic ex- 


pansion, showing that it represents a diverging wave. In this expansion 

the factor (pm/c) r+ ^ in the denominator yields (W|) in the 

interpretation, and the extra II function in the denominator often saves 
the convergence. 


8.51- A wave spreading out symmetrically from a point source gives 
a velocity potential depending on the operator 



Now suppose that a disturbance occurs at time 0 all along the axis of 
z, giving a disturbance of the form (1) per unit length. Then if f is the 
^-coordinate of the place where the disturbance starts, the disturbance 
at a point on the plane z = 0 distant m from the axis is 



<£- f - e~ prla dC 

j -00 r 

(2) 

Now put 

r — 'us cosh v, £ = m sinh v 

(3) 

Then 

<j> = 2 j exp cosh dv 



«*{?)•- 

W 


If the intensity of the source is proportional to cos K we write on the 
plane z = 0, 

4>= ( * cos 

J- 00 V 

-2 J cos (A.zet sinli^) exp c°sh (5) 

Split the cosine into exponential terms and put 

= y COS a; Xw = y sin a. (6) 


♦-jCV 


cobIi (t>-ca) 4. coflli {«+ia)j. 


Then 
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which is independent of a for — Jir < a < and therefore 
= 2 J e~ ycmhv dm 

= ^ 0 {(| + X»)*«r} (8) 

For a point with coordinate z we can put 

(9) 

and cos A£ = cos Xz cos Xz' — sin Xz sin Xz' (10) 

The point is on the plane z = 0, and we can repeat the process; evidently 
by symmetry the term in sin Xz f contributes nothing, and 

<£ = 7 r JTo + A. 2 ^ gt| cob Xz (11) 

If the source has unit intensity between £ = ± h and zero elsewhere, 
we can write the intensity as 

± f {eMf+M-eMf = e^sinhXA — 

2ttlJ l 1 X mj L X 

2 r . .jdX . , 

= - / cosA£smAA — . ...(12) 

Thus 0 = 2 jT cos Xz sin Xh K n |(^ + \ 2 ^ m j ^ (13) 

If h becomes very small we have for the disturbance due to unit intensity 
over the short range + h f 

= 2k J cos Xz jfir 0 + A. 5 ^ t5t| dX (14) 

8.52. We can put 8.51 (1) in an alternative form in cylindrical 
coordinates. On the axis of z , with 2 ? positive, we have 

I 0-* r /° = - e-« s /o = f 0 -^dK (1) 

9* Z JkIc 

Also the solution of the equation 

v a <£ = '^r (2) 

that is equal to e~ Ka when wr - 0 is e~ Ka J 0 j^A 3 - j. Hence 

l e ~ Kr ' c= [ lc e ^ J ^{( xl -^ w } dK 

= - kdk , 

Jo (W+f/c*)* 

whence ^ (f - ^ is expressible in cylindrical coordinates. 
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8.6. The Submarine Gable . The submarine cable is a uniform con- 
ductor with self-induction, capacity, resistance, and leakage. At distance 
sc from the end suppose that the charge per unit length is y, the potential 
V, and the current i. Then the following equations hold : 




l dt + “ 


dV 

das’ 


( 1 ) 

( 2 ) 


dy 

d t 


sV. 


( 3 ) 


l, k, r are the self-induction, capacity, and resistance per unit length. 
The leakage is such that a potential V produces a current sV per unit 
length escaping through the insulation. Up to t = 0, y, i, and V are zero; 
afterwards V is raised to V oi which may be a function of i £, at # = 0. 


Then the subsidiary equations are 

(l p + r)i = jr. (kp + s) V=- S £., ( 4 ) 

d 2 V 

whence ^ = (lp + r) (kp + s ) V. (5) 

Put (Ip + r) ( kp + s) = q 2 = Ik {(p + p) 2 - cr 2 } (6) 

Then the operational solution is 

F= e -<i‘F 0 (7) 


If the self-induction and leakage are negligible we can put £=0, s = 0, 
and q 2 - krp . Then the solution has the same form as for conduction of 
heat. This condition occurs in ordinary telegraph wires. If in this case 

K=im 


V= 1 - erf 


.( 8 ) 


2 (krift 

In the general case we put 

Ik- l/c a ..(9) 

Since p + <r and p - o- are equal to r/l and s/k, which are both positive, 
p is positive. We can without ambiguity take or positive; then c r<p. 
If one of r/l and s/k is zero, o- = p; if both are equal, <r = 0. These are 
the extreme cases. 

If cr = 0 we have simply 


r=exp{-(^ + P )?}.F„ 


= <rp»yo y. 


.( 10 ) 


BESSEL EXTINCTIONS 


105 


Then the variation of V with time at distance x from the origin is an 
exact copy of that at the end except for the delay cc/c and the constant 
attenuation factor e~~v xtc . This is the distortionless cable . 

If a is not zero we may take 

Vo = H(t), 

and the solution is 


F=exp[-|{( i? +p) 2 -^}4] ( 11 ) 

(«) 


The form (11) can be treated by expanding in powers of cr 2 ; the solution is 






.( 13 ) 


the terms, which are all positive, rising in powers of (^“f) an( ^ 

o- 2 (t - . So long as these numbers are small, this series may be use- 

ful. But if o-xjc is large the series is not a useful form for any appreci- 
able range of t after the arrival of the disturbance. 

The form (12) may be treated by putting 

^ = ( 14 ) 


Then f =£* (15) 

The integrand has branch points at /x = ±<r and a pole at /x = p. If we 
write 


/ 0) (a* 2 - O' 2 ) 4 . (!6) 


/'(/*) = *- 
/"60 = 


JlX 


C - <J-f ’ 

arx 


C (/x 3 — cr 2 )^ 

then f (/x) vanishes at /x = /x 0 , where 

•"( i -&)"*• 


/x 0 = 


( 17 ) 

.( 18 ) 


( 19 ) 


Whenever t > x/c, /x 0 is greater than <r. If t is just greater than xjc , /x 0 
is large \ if t is large, /x 0 approaches cr, which is in general less than p. 
Hence for values of t greater than some critical value /x 0 is less than p, 
and a path through p 0 parallel to the imaginary axis will have p on its 
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right, and. the pole there must be allowed for separately. This value is 


given by 

ctjx = (1 - o s /p‘)~^ (20) 

We have /(^>) = <rt (l - , (21) 

™-3(: i-£)'. W 


On a path through p 0 parallel to tlie imaginary axis, considerable values 
of the exponential are spread over a range of order {/"OOF^ or 

(bO* (* - ??) * b8 " B to * tbe ('-»)*• 

which is large if t is just greater than oojc but becomes small if t is large. 
In the former case we cannot neglect the variation of p. - p on the path 
of steepest descent. In the latter, the ratio is much less than (o-#/c)“^ 
which we have already assumed small. Then the contribution from the 
neighbourhood of p 0 is 


2x(p-^) 0Xp/W -{/" 2 (^)} 


( X (i 

a,- 2 V 

\2ircW V 

any 


-f r 




x exp {<rt (l-~) -pt} 


.(23) 


But the pole at p is to the right of the path and makes a further 
contribution 

0 - 0 * exp jp# - ® (p 2 - cr 2 )^ = exp j- ® (p 2 - cr 2 ) J |. . . .(24) 

This is independent of the time. The expression (23) tends to zero when 
t is large, since p>cr. The sum of (24) and (23) is an asymptotic form 
of the solution when t is greater than the value given by (20) ; and the 
potential approximates to the steady value (24) if ( ct/$J exceeds 
(1 - o- 2 /p 2 ) -1 . Thus the time needed for the disappearance of the variable 

part is of order ~ (l - 

In the critical case where either the resistance or the leakage is zero, 
<r = p, and this time is infinite, so that the value (24) is never approached. 
We can obtain a solution in a series from (15) by writing 

p“ — o’ 8 = (2v — /a) 8 . 


< 25 ) 
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Then P = v + £i 2v ~P = ''-£ v ( 26 ) 

e -p * r f t x \ v 2 / a?\i 4 v*-o* dv 
V= K iz eXP 1 V ~ o) ' + 4v (* + ;)} 4 ^4^ 7 

(27) 

exp ( vtl + i t *)( 1+ 7=7 1 + T * (28) 

where tfi = £ 2 = tf + -) (29) 

C 0 

2vi = p + (p a - a- 2 )^; 2v a = p — (p 2 — cr 2 )^ (30) 

But we have the general expansion 

e*p{j»(**s)}-J_/-4M, (81) 

whence 

-»(* + 5*h .1 ££-2) , }‘* {'(<■-?)*}■ <*» 

The integral in (28), when t x is positive, may be taken around any circle 
large enough to include the singularities at v = 0, v u v 2 . But negative 
values of n make no contribution to the contour integral ; whence 

(33) 

But ( v n — = 0 if m>o) , 

and =1 if n - 0 J 

> U .*.* > fci 

27TlJc V — I'l v 27Tt J r= l \V / V >. (ooj 

and = 0 if « = 0 J 

Hence 

r= e -» < [i;,{<r(« l # a ) 4 }+J i {|@ i f«+^)4{^(#,# a ) i }]- (36) 

In the case where cr = 0, each of the (2 /o , ) n / n {<r(£i£ a )fy reduces] to its 
first term {(Ma)^} 71 /^ Vi = p, v 9 = 0, and 

(37) 

agreeing with (10). 
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Alternatively if we replace 3/3#! operationally byjpu (28) is equivalent 
to 


nl o*t t 

e ' piex P^ 


■( 


1 + 2 
r=l 


vj+vj 

Pi 


)h { o 




- J 2 1 f^Y +22 W < + dy, >, 

~ L=0»! V4ft ) n=0 r=l » 1 4” ft** / " ^ 

_- P< f ; (WO*. , w 5 QgW if/ y f , ,.,0 

U=0 w! »! »=o r=i n\ (n + r)\ 1 '* *ty’ 

which on rearrangement gives the form (36). 


(38) 


8.7. The operational method raises numerous points concerning 
the solution of the differential equations of mathematical physics. 
Chapter I of this work suffices to show that it gives the correct results 
when the system has a finite number of degrees of freedom ; but as soon 
as we proceed to continuous systems we find that the solution involves 
operators not expressible by series of negative powers of p, and the 
method of Chapter I is no longer available as a justification. 

We have seen that in certain continuous systems the use of the partial 
differential equation is itself open to suspicion. Thus in a string the 
value of dyjdx may be discontinuous at certain points, and dh/fdx 2 does 
not exist there and the differential equation is simply meaningless. 
The real problem that concerns us is not whether the operational method 
gives the correct solution of the partial differential equations, for those 
equations are themselves defective. The question is whether, if we use 
these equations, modified to allow for the initial conditions, and solve 
them by the operational method, the answer is correct for the actual 
physical system. We have to test the answer in relation to the physical 
problem, not in relation to the differential equations, which are them- 
selves an imperfect statement of the problem and in any case only an 
intermediate stage in its solution. My own view is that most continuous 
systems should be regarded as the limits of discrete systems. Thus if 
we have a string loaded at intervals l with particles of mass pl } and the 
displacement of the rth particle from the end is y r> the equation of 
motion of a particle is 

Hr = ~ ja (2jfr ~ j/r-l ~ 2#r+l) > 
where <f=P/p, P being the tension. 


( 1 ) 
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If we make l approach, zero while p and c remain constant, the system 
approaches in the limit the continuous string with line-density p. If 
further the quantity on the right tends to a definite limit we can put 
rl - x, so that x is the distance of the particle from the end, and the 
equation takes the ordinary form 


dt 2 dtf • 


( 2 ) 


But whether the limit 9 2 y/9# a exists or not, (1) represents the equations 
of motion of a system with a finite number of degrees of freedom, and 
for any initial conditions there is an operational solution involving only 
powers of p~\ and therefore certainly equivalent to the corresponding 
integral along the line L. It actually involves operators of the form * 



which can clearly be expanded in negative powers of p . But the 
corresponding function of X, when l becomes very small, is for large 
values of r 


<j> (\) = exp 
= exp 


[-2rlog{l + y+o(^)}] 

[-{b^S 


and therefore tends to e~ rKl/o or Thus the limit of the solution is 
precisely Bromwich's interpretation, with (p) replaced by its formal 
limit e~ pxl °. The operational solution is therefore justified as the limit 
of the actual solution for a discrete system. 


Jeffreys, Proc. Gamb. Phil Soc. 23, 1927, 768-778. 


NOTE 


ON THE NOTATION POE THE ERROE FUNCTION 
OR PROBABILITY INTEGRAL 

The notation Erf x is one that I adopted in 1916 under the impression 
that it was in general use, and I haye since used it in several publications*. 
My definition was queried by a correspondent in 1925, and I have not suc- 
ceeded in tracing its origin. I have, on the other hand, discovered a surprising 
confusion of other notations. The earliest, due to Gauss t, is 

e (*)— L [\-*dz. (i) 

>JirJ 0 

No name is given to the function by Gauss, and there is no sign that he meant 
the notation to be permanent. Of modern writers, Carslaw, Brunt, and Coolidge 
use this notation. Fourier { gives 

-L r e-r % dr, (2) 

sjir J R 

but has, so far as I have traced, no modern followers. Jahnke and Emde, in 
their tables, use 

$(#)= ~= [ e - * 2 dx, (3) 

sJtt J o 

and call this the Fehlerintegral. The same notation is widely used by other 
German writers. Whittaker and Robinson § also use it and call the function 
the Error Function. The notation 



was introduced by J. W. L. G3aisher||, who also used IT 

Erfc x = j e~ dx—^Jir — Erf x (5) 

The latter function is also called erf x by R. Pendlebury** 

* Phil, Mag. 32, 1916, 579-585; 35, 1918, 273; 38, 1919, 718. M.N.B.A.8. 77, 
1916, 95-97. Proc . Boy. Soc. A, 100, 1921, 125-6. The Earth , 1924, 1929, 
t Werke , 4, 9. First published 1821. 

X TMorie Analytique de la Ghaleur , 1822, 468. 

§ Calculus of Observations , 1924, 179. 

Ji Phil. Mag . (4), 42, 1871, 294-302. 

IT Loc. cit . 421-436. 

** Loc . cit. 437-440. 
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Whittaker and Watson* use Glaisher’s notation with the mannings of Erf 
and Erfc interchanged, while Jeans’st erf x is — times Glaisher’s Erf. 

sjn 

Numerous other writers use the integrals, but omit to give any special symbol 
or use only the non-committal “I.” Heaviside + used erf a to denote the 
function (1) or (3). 

It can hardly be denied, in view of the wide application of the error function 
to thermal conduction, the theory of errors, statistics, and the dynamical theory 
of gases, that it merits a distinctive notation. It is equally clear that none of 
those yet used has obtained general acceptance. Of them, it seems that those 
involving only the single letters 6, * are wholly undesirable. In dynamical 

problems these letters are in continual use for couples, while <£ is often wanted 
for a velocity potential and ^ or f for a stream function ; ^ is of course also 
often needed for an angle. Further, © has an established meaning already in 
the theory of elliptic functions, namely, the Theta Function of Jacobi, while 
<I> has another meaning in the theory of numbers. To avoid confusion in some 
of the applications of the function it seems necessary to use a combination of 
letters, and the only one in frequent use is Erf. On the other hand, as 
Dr Goldstein has pointed out to me, the capital letter is somewhat incon- 
venient in typing, and is not intrinsically useful, so erf has some advantage. 

As to what function should be denoted by this symbol, the most convenient 
is indicated by the practice of the compilers of tables, such as Jahnke and 

Emde and Dale, who tabulate 


=/■ 
7T J 0 


e~ u du. This is an odd function and 


becomes unity when #=co , two properties that make for analytical convenience 

and are connected with the fact that this form, including the factor % r ” i, 
usually occurs as such in the solutions of the relevant problems. Accordingly 
I think that, of the various notations proposed, 


erf.- ■L/* 

JttJo 


e~ u du 


is the most convenient, and is worthy of general adoption. It has also the 
recommendation that it has previously been used by Heaviside. 


* Modern Analysis, 1915, 335. 
t Dynamical Theory of Gases , 1921, 84. 
X Electromagnetic Theory , 2, 51. 


INTERPRETATIONS OF THE PRINCIPAL OPERATORS 


p~ n H(t )= 0 t< 0 
! 


<> 0 . 


p-a ’ ( p — a) n (?i — 1)1 


,<*£ , 


-2 — e a< — 1. 

p — a 


Jp=^+-/~^ anrt 
f(p) m , y m e at 
F(P) F(0) + faF‘ia) e ‘ 

Jr~i= sin nt > -.,~ a =oos«i; -^-5=amh?i<i 
p‘+n- p i +7i i ’ p* - ri* 


p*-w 


,=cosh nt. 




V 


li p— h 2 q'\ 

e~1*= 1 - orf — ; 2 e-«*=— e ~^ im ; 


- «-«' r =2A (iV o~ x ' wt -x (l - orf -£-} ; 

? W V 2/t<V 

— =l-arf -L_oxp (aVA+aa.’) -fl-orf +a/rfi')l. 
-« a/rfl l '2 At* J> 


ae 


P 


{p + (p 2 - l)fy* (p 2 - 1 


=• 4 ( 0 ; 


P 


{p+(p 2 + i) 4 } n (p 3 +i) 4 


= 4 ( 0 - 


where 


TABLE of e~ x \ erf x, 1 - erf x, <E>! (a?), <2> a (a?), 


<E> X (a?) = / (1 — erf a?) da?, $2 (a?) = f ^1 (a?) da?. 
J X J X 


£ 

«-* 

erf x 

1 - erf x 

*|M 


o-o 

1*0000 

0*0000 

1*0000 

0*5641 

0*2500 

0*1 

0*9990 

0*1125 

0*8875 

0*4698 

0*1984 

0-2 

0*9608 

0*2227 

0*7773 

0*3866 

0*1556 

0-3 

0*9139 

0*3286 

0*6714 

0*3142 

0*1206 

0-4 

0*8521 

0*4284 

0*5716 

0*2521 

0*09241 

0*5 

0*7788 

0*5205 

0*4795 

0*1996 

0*06990 

0-6 

0*6977 

0*6039 

0*3961 

0*1559 

0*05220 

0*7 

0*6126 

0*6778 

0*3222 

0*1201 

0*03846 

0*8 

0*5273 

0*7421 

0*2579 

0*0911 

0*02795 

0*9 

0*4449 

0*7969 

0*2031 

0*0682 

0*02003 

1*0 

0*3679 

0*8427 

0*1573 

0*0502 

0*01415 

1*1 

0*2982 

0*8802 

0*1198 

0*0364 

0*00985 

1*2 

0*2369 

0*9103 

0*0897 

0*0260 

0*00675 

1*3 

0*1845 

0*9340 

0*0660 

0*0183 

0*00456 

1*4 

0*1409 

0*9523 

0*0477 

0*0126 

0*00302 

1*5 

0*1054 

0*9661 

0*0339 

0*0086 

0*00197 

1*6 

0*0773 

0*9763 

0*0237 

0*0057 

0*00126 

1*7 

0*0556 

0*9838 

0*0162 

0*0038 

0*00079 

1*8 

0*0392 

0*9891 

0*0109 

0*0024 

0*00049 

1*9 

0*0270 

0*9928 

0*0072 

0*0015 

0-00029 

2*0 

0*0183 

0*9953 

0*0047 

0*0009 

0-00017 

2-1 

0*0121 

0*9970 

0*0030 

0*0006 

o-oooio 

2*2 

0*0079 

0*9981 

0*0019 

0*0003 

0-00005 

2*3 

0*0050 

0*9989 

0*0011 

0*0002 

0*00003 

2*4 

0*0031 

0*9993 

0*0007 

0*0001 

0*00001 

2*5 

0*0020 

0*9996 

0*0004 

0*0000 

0*00001 

2*6 

0*0012 

0*9998 

0*0002 


0*00000 

2*7 

0*0007 

0*9999 

0*0001 



2*8 

0*0004 

0*9999 

0*0001 



2*9 

0*0002 

1*0000 

0*0000 



3*0 

0*0001 

1*0000 

0*0000 
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